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1 Introduction 



1.1 Goals of W-pair production 

The measurements performed at LEPl have provided us with an extremely accurate knowledge 
of the parameters of the Z gauge boson: its mass, partial widths, and total width. Except per- 
haps for the bb and cc decay widths all data are in perfect agreement with the Standard- Model 
(SM) predictions 0, |]. There even is first evidence that the contributions of gauge-boson 
loops to the gauge-boson self-energies are indeed required [Q. Thus an indirect confirmation of 
the existence of the triple gauge-boson couplings (TGC's) has been obtained. 

The future measurements of W-pair production at LEP2 will add two important pieces of 
information to our knowledge of the SM. One is a determination of the W mass, which at 
present is only directly measured at hadron colliders. The envisaged precision of 40-50 MeV 
gives a significant improvement on the present Tevatron measurements Since the W mass 
is one of the key parameters of the electroweak theory, such an improved accuracy makes the 
tests on the SM more stringent. 

The second piece of information that W-pair production can provide is the structure of the 
triple gauge-boson couplings. These couplings now play a role in the tree-level cross-section, in 
contrast to LEPl physics where they only enter through loop corrections. So the Yang-Mills 
character of the TGC's can be established by studying W-pair differential cross-sections at 
LEP2 energies. Since these couplings are at the heart of the non-abelian gauge theories, this 
information is essential for a direct confirmation of the SM. 



1.1.1 Measurement of the mass of the W boson 

With the precise measurement of the mass of the W boson at LEP2 the situation for the 
electroweak input parameters changes with respect to LEPl. The common practice at LEPl 
is to use for Mw a value derived from the Fermi constant G^, which is accurately known 
from muon decay. The relation to obtain Mw follows from the SM prediction for muon decay 



' V2M^{1-M^/Mi)l-Ar' 

where Ar = at tree level and where Ar is rrit- and Mn-dependent when loop corrections 
are included. Thus, Mw in LEPl calculations is rrit- and Mn-dependent through the above 
procedure. At LEP2, where one wants to measure Mw and, hence, wants to treat Mw as a 
fit parameter, the above relation now primarily acts as a test of the SM. The above relation 
predicts for any chosen Mh and measured Mw a value for rrit that can be used as input for 
LEP2 loop calculations and can also be compared with the directly observed top-quark mass 
from the Tevatron 0. 
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As to the actual procedure to measure M-w from W-pair production, two methods are 
advocated |^. A third one, involving the measurement of the endpoint energy of a charged 
lepton originating from a decaying W boson, was suggested a few years ago but turns out to 
be less promising than other methods. 

One procedure requires a measurement of the total W-pair cross-section close to threshold, 
where the size of atot is most sensitive to the W mass. The energy proposed is 161 GeV. For 
this method the theory should obviously predict Utot with a sufficient accuracy, to wit ~ 2%.Q 
Therefore the radiative corrections (RC's) to the total cross-section should be under control. 

The other method looks at the decay products of the W bosons, in particular at the four- 
and two-jet production. From the measured momenta of the decay products one tries to 
reconstruct the W mass. In this reconstruction a good knowledge of the W-pair centre-of- 
mass energy is essential. Since there will be an energy loss due to initial-state radiation (ISR) 
of primarily photons, the W-pair energy will be different from the laboratory energy of the 
incoming electrons and positrons. So, for this method the ISR should be well under control, 
i.e. the theoretical error on the average energy loss, {E^) = J dE^ (da/dE^) E^/atot, should 
be less than 30MeV, which translates into a theoretical error of less than 15MeV on the 
reconstructed W mass. This again is an aspect of radiative corrections. 



1.1.2 Test of non-abelian couplings 

Within the SM the triple gauge-boson couplings have the Yang-Mills (YM) form. Amongst 
others, this specific form for the TGC's leads to a proper high-energy behaviour of the W-pair 
cross-section and is a requirement for having a renormalizable theory. Couplings different from 
the YM form, called anomalous or non-standard couplings, will in general lead to a high-energy 
behaviour of cross-sections increasing with energy and thereby violating unitarity. At LEP2 
the energy is too low to see such effects and in order to establish the presence of anomalous 
couplings one therefore has to study in detail the angular distributions of the W bosons and their 
decay products. In particular, the angular distribution in the W-production angle 9 is sensitive 
to non-standard couplings. Again, the knowledge of RC's to the tree level SM predictions is 
required, since they affect the Born-level angular distributions 



As elsewhere in this volume a detailed report on non-standard TGC's is given |T^, only a 
few comments will be made here. 

When one considers the most general coupling between three gauge bosons allowed for by 
Lorentz invariance, assuming the gauge bosons to be coupled to conserved currents, one ends 
up with seven possible couplings for the ZWW and 7WW interaction. Of these seven there are 



three which are CP violating and one which is CP conserving but C and P violating [11 



"Throughout this report the reqmred theoretical accuracy is taken to be half the expected statistical error. 
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In practice it will be impossible to set limits on all these couplings. Therefore usually some 
assumptions are made to reduce the number of 14 couplings [|ll|, |12|. For instance, one may 
restrict oneself to CP-conserving non-standard couplings so that 8 couplings are left. Of these 
the electromagnetic ones can be reduced further by omitting the C and P violating one and 
requiring the strength of the electromagnetic coupling to be determined by the charge. Two 
possible anomalous electromagnetic couplings remain and four ZWW anomalous couplings. 
Even with this reduced number it will be impossible to set experimental limits on all of them 
simultaneously. 

However, there are theoretical arguments that such a purely phenomenological approach 
is also not required. First of all one might use symmetry arguments, motivated by specific 
models for the non-standard physics, to find relations between the anomalous couplings . 



Alternatively, when one considers the electroweak theory as an effective theory originating from 
a field theory that manifests itself at higher energies, then also some small anomalous couplings 
may be present at lower energies. In such a SU{2) x U{1) gauge-invariant framework the non- 
standard physics, situated at an energy scale A, decouples at low energies and the anomalous 
TGC's are suppressed by factors {E/AY~'^, according to the dimension (d) of the corresponding 
operators. This naturally introduces a hierarchy amongst the anomalous TGC's based on the 
dimension of the corresponding operators [p^ . 

From the perspective of this report the origin of the non-standard couplings is not so im- 
portant, but the fact that they often modify angular distributions is relevant. Also SM effects 
- like RC's, the finite decay width of the W bosons, and background contributions - provide 



deviations from the tree-level distributions [Bl 10 . 



1.2 How to obtain accurate predictions? 

In order to extract the wanted information from W-pair production it is clear from the remarks 
above that RC's are needed for total and differential cross-sections and, moreover, for the de- 
termination of the energy loss. Anticipating ~ 10^ W-pair events, the theoretical accuracy 
that should be targeted for the SM predictions is ~ 0.5%, although specific final states, distri- 
butions, or observables in fact often do not require such a precision (like, e.g., atot at 161 GeV 
or the energy loss). 

Ideally one would like to have the full RC's to the final state of four fermions, which originate 
from the two decaying vector bosons. In practice such a very involved calculation does not exist 
and is hopefully not required in its complete form for the present accuracy. For the discussion 
of the LEP2 situation and strategy it is useful to distinguish three levels of sophistication in 
the description of W-pair production. 

The first level is to consider on-shell W-pair production, e^e~ — >■ W"^W~, which at tree 
level is described by three diagrams: neutrino exchange in the t channel, and 7 and Z exchange 
in the s channel. Here the complete 0{a) RC's are known, comprising the virtual one-loop 
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corrections and the real-photon bremsstrahlung [n, 15 1. When one wants to divide the 0{a) 



corrections into different parts the situation differs from LEPl [^. 

As to the bremsstrahlung, a gauge- invariant separation into initial-state radiation, final-state 
radiation (FSR), and its interference is not possible like, e.g., in /i-pair production at LEPl. The 
reason is that the photon should couple to all charged particles in a line of the Feynman diagram. 
The t-channel diagram then makes a separation into ISR and FSR impossible. However, the 
leading logarithmic (LL) part of ISR is in itself gauge-invariant. This can be combined with 
the LL QED virtual corrections so that a LL description with structure functions for ISR can 
be given |l^ . 



A separation of the virtual corrections into a photonic and weak part is also not possible 
since charged vector bosons are already present at Born level, necessitating an interplay between 
7- and Z-exchange diagrams in order to preserve SU{2) gauge invariance. 

Once one has a description of on-shell W-pair production one can attach to it the on-shell 
W decay. Again, the RC's to this decay are known ||T7[|-||20[|. 



The next level of sophistication is to consider off-shell production of W pairs, which then 
decay into four fermions ||21|. The W propagators with energy- dependent widths can be taken 
into account. Although this description of four-fermion production through virtual W bosons 
is a natural extension of the on-shell evaluation, it is not a gauge- invariant treatment. In fact 
there are more diagrams needed to calculate such a four-fermion process. As to RC's, the ISR 
and FSR can be implemented in LL approximation, but the full set of virtual corrections have 
not yet been calculated. 

The final level for the study of W-pair production would be a full 0(Q;)-corrected evaluation 
of all possible four-fermion final states. At tree level there now exist evaluations where, besides 
the three off-shell W-pair diagrams, all other diagrams for a specific four-fermion final state 
have been included p2|. On top of that ISR can be taken into account. Again, one has to 



be aware of possible gauge-invariance problems. In particular, the introduction of energy- 
dependent widths in the W-propagators will destroy electromagnetic gauge invariance and may 
introduce dramatically wrong cross-sections in certain regions of phase space. 

As long as a full (9(Q;)-corrected evaluation of four-fermion production is not available, 
certain approximative schemes, like for instance the 'pole scheme' [^-[^, may be useful. This 



goes beyond the treatments where only ISR through structure functions is taken into account. 
Actual numerical results from a complete 'pole-scheme' evaluation are not yet available. 



2 On-Shell W-Pair Production and W Decay 

Although the actual process that will be probed at LEP2 is e+e~ —>■ 4f['j,g), we first focus on 
the production and subsequent decay of on-shell W bosons, being basic building blocks in some 
of the schemes for handling off-shell W bosons. In contrast to off-shell W-pair production the 
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on-shell processes are not plagued by the problem of gauge invariance for unstable particles and 
a complete set of 0{a) radiative corrections is available. Consequently, they are well suited for 
studying the typical sizes of various radiative corrections. Moreover, many of the important 
features of the production and decay of off-shell W bosons are already contained in the on-shell 
limit. We indicate explicitly where the width of the W bosons radically changes the on-shell 
predictions. 



2.1 Notation and conventions 

We use the Bj0rken-Drell metric g'^jy = diag(-|-l, —1, —1, —1) and fix the totally antisymmetric 
tensor by e^-*^^^ = The matrix 75 is defined as 75 = i7°7^7^7^ and the helicity projectors 
u!±, which are used to project on right- and left-handed massless fermions, as 

a;± = ^(l±75). (2) 

First we set the conventions for the process 

e+(p+, + e-(p_, «_) ^ W+(A;+, A+) + W-(A;_, A_), (3) 

where the arguments indicate the momenta and helicities of the incoming fermions and outgoing 
bosons {Ki — ±|, Aj = 1, 0, —1). Note that we sometimes use the shorthand version = ± in 
certain sub- and superscripts. In the centre-of-mass (CM) system of the e+e~ pair, which we 
will refer to as the laboratory (LAB) system in the following, the momenta read 

= E(l, 0, 0, Tl), = E(l, T/3sme, 0, T/?cos^), (4) 

with E denoting the beam energy, the scattering angle between the e"*" and the W"*", and 



[3=^1 — M^/E"^ the velocity of the W bosons. The Mandelstam variables used in the 
following are given by 

s = (p+ +p_f = {k+ + k^f = ^E\ 

t =[p^-k+f = {p_-k_f = -E'^ (I + [3^-2(5 cofi 6), (5) 
u= {p+-k-f = {p-- k+y = -E\l + (3^ + 2(5 cos 6). 

In order to define helicity amplitudes we need to introduce the corresponding polarization 
vectors for the W+ and W~ boson'' 

£±^(^±,+1) = ^(0, =Fcos^, -i, ±sin6'), 
v2 



^Note that the hehcity of the massive W bosons is not Lorentz-invariant. We define it in the LAB system. 
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^±''(^±,-1) = ^(0, =Fcos6', +z, ±sin6'), 
V 2 

£±^'(^±, 0) = -^(/3,Tsin^,0,Tcos^). (6) 

Because we are neglecting the electron mass, the helicity of the positron is opposite to the 
helicity of the electron 

K_ = — = K. (7) 

Henceforth we refer to the helicity amplitudes for W-pair production as 7V1 (k, A+, A_, s, t). CP 
invariance implies the relation: 

M{k, A+, A_, s,t) = M{k, -A_, -A+, s,t). (8) 

This holds in the SM if we neglect the CP-violating phase in the quark-mixing matrix. Even 
in the presence of this CP-violating phase, the CP breaking occurs at 0(q!^) in SM W-pair 
production and is additionally suppressed by the smallness of the mixing angles between the 
quarks. Consequently in the SM one effectively has only 12 independent helicity-matrix ele- 
ments instead of 36. In the presence of substantial (non-standard) CP violation this number 
should be increased to 18. This allows for a decomposition of the matrix elements in terms of 
an explicit set of 12 (18) independent basic matrix elements multiplied by purely kinematical 
invariant functions (coefficients) |14|, |1^. 

From the helicity amplitudes the differential cross- sect ions for explicit W-boson polarization 
and various degrees of initial-state polarization can be constructed. For example the differential 
cross-section for unpolarized electrons, positrons, and W bosons is given by 

For the numerical evaluations we have to fix the input parameters. We use the following 
default set |,|2|]: 



a = a(0) = 1/137.0359895, = 1.16639 x 10"^ GeV"^ 

Mz = 91.1884 GeV, Mr = 300 GeV, 

me = 0.51099906 MeV, = 105.658389 MeV, = 1.7771 GeV, 

mu=47.0MeV, mc= 1.55 GeV, 

md = 47.0 MeV, = 150 MeV, = 4.7 GeV. 

The masses of the light quarks are adjusted in such a way that the experimentally measured 



hadronic vacuum polarization is reproduced. In the actual calculations either these light 



quark masses are used or the dispersion-integral result for the hadronic vacuum polarization 
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Mh [GeV] = 


60 


300 


1000 


Mw [GeV] 


mt [GeV] 


80.10 


119.1 


137.3 


154.2 


80.18 


133.9 


151.6 


168.0 


80.26 


148.1 


165.3 


181.2 


80.34 


161.7 


178.3 


193.9 


80.42 


174.4 


190.7 


206.1 



Table 1: Calculated mt for a;s(M|) = 0.123 and different Higgs- and W-boson masses, using 



the state-of-the-art calculation described in section |2.3.3| . The theoretical error in mt is roughly 
1-2 GeV. 



of |27|. The strong coupling constant is calculated according to the parametrization of | 28| , 
using as(M|) = 0.123 as input. The W-pair threshold region is very sensitive to and 
consequently, if Mw were to be calculated from the other parameters using the muon decay 
width including radiative corrections (like at LEPl), to the masses of the top quark and Higgs 
boson. This is of course not very natural, since we want to use Mw as a model-independent 
fit parameter. In view of this we use a, G^, Mz as input^, treat Mw as free fit parameter, and 
calculate mt from muon decayQ. This calculated value of mt can then be confronted with the 
direct bounds from Fermilab (weighted average mt = 180 ± 12 GeV [Q) and the indirect ones 
from the precision measurements at LEPl/SLC (mt = 180^9^20 GeV 0). In this scheme mt 
hence plays the role played by Mw at LEPl. The above set of parameters and the default W 
mass Mw=80.26 GeV g, |] yield 

mt = 165.26 GeV. 

At present the error on this value as a result of the 160 MeV error on Mw is roughly 27 GeV. 
The precise Mw measurement at LEP2 will lead to a reduction of this error by a factor of 
four. Of course the so-obtained top-quark mass will become as- and Mn-dependent through 
the RC's associated with the muon decay. We come back to that point in section |2.3.3 and 



only mention here that the Mh dependence of mt amounts to roughly 33 GeV in the range 
60 GeV < Mh < 1 TeV, which clearly exceeds the above-mentioned expected error from Mw 
(see table |1|). 

Finally, the sine and cosine of the weak mixing angle are defined by 

Cw = cos6'w = sw = sin6'w = \/l - c%j. (10) 



'^We do not eliminate Mz or G^, as Mz is needed in the RC's and reduces the size of the RC's to the 
production and decay of the W bosons. 

''it should be noted that the RC's associated with the muon decay make that this parameter set is not 
overcomplete. 
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Figure 1: Lowest-order diagrams for e+e W+W . 
2.2 Lowest order 

In the SM there are three lowest-order diagrams (figure |l|), if we omit a Higgs-exchange diagram 
that is suppressed by a factor rrte/Mw and thus completely negligible. The t-channel diagram 
involving the exchange only contributes for left-handed electrons. The s-channel diagrams, 
containing the non-abelian triple gauge-boson couplings, contribute for both helicities of the 
electron. The corresponding matrix element reads 



e2 1 

MboUh, A+, A„, s, t) = - Ml + 



s sw s - Ml 



2{Ml-M^^), (11) 



with = 1 for left-handed electrons and 5^- = for right-handed electrons, and 

M'l = v{p+)^+{k+,X+) (|^+ ^_(fc_,A_)^«M(p_), 

= v{p+) t±_JL_ \e^{k^^ A+) ■ £_(A;_, A_)] ^.m(p_). 

Ml = v{p^) {UiK. A+) [^4^-, A_) ■ - ^_(fc-, A^) [£+(A;+, A+) ■ uj^P^)- (12) 
After inserting the explicit form of the Z-boson-fermion couplings 

« "^W 1 /-ION 

9eeZ = ^^-7^ > (13) 

we can organize the lowest-order amplitude into two gauge-invariant subsets: 

A4Bom(«:, A+, A_, s, t) = Mi{k, A+, A_, s, t) + e^Mgin, A+, A_, s, t), (14) 

where 

MiiK, A+, A_, s, t) = i + 2 (A^^ - A^^), 



Mq{k, A+, A_,s,t) 



1 



s s - M| 



2{Ml-M^,). (15) 



The gauge invariance of the two contributions Mi and Mq can be simply inferred from 
the fact that they are accompanied by different coupling constants, one of which involving the 
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electromagnetic coupling constant e, the other the charged- current coupling constant e/(v^sw)- 
Whereas Adj corresponds to the pure SU{2) contribution, the parity-conserving contribution 
A^Q is a result of the symmetry-breaking mechanism. 

The lowest-order cross-section determines the essential features of W-pair production. The 
threshold behaviour is important for the determination of the W mass from the measurement 
of the cross-section in a single energy point very close to threshold i.e. at -y/i = 161 GeV. 
For small 13 the matrix elements behave as 

A^^,A^^oc/5, A^^ocl (16) 

for fixed scattering angles. Consequently, the s-channel matrix elements vanish at threshold^ 
and the t-channel graph dominates in the threshold region. For /5 -C 1 the differential cross- 
section for unpolarized beams and W bosons is given by P] 




1 . .3c^-l 

Ad 



Born ^ ^-^W 



l + A^cosO^ + 0{I3 



(17) 



where the leading term oc (3 originates from the t-channel diagram only. Note that the leading 
term is angular- independent. In the total cross-section 

aBorn^^TV4/3 + 0(/3=^), (18) 
S 

all terms oc jS"^ drop out and the s-channel and the s-t-interference contributions are propor- 
tional to (3^. This is the consequence of CP conservation, fermion-helicity conservation in the 
initial state, and the orthogonality of different partial waves 0. Hence in the threshold region 
the t channel is dominant and the cross-section for e+e" — >■ W+W" is not very sensitive to the 
triple gauge-boson couplings. 

In table ^ we give the integrated cross-section for different centre-of-mass energies and dif- 
ferent polarizations of the electrons (+,—). Positrons are assumed to be unpolarized. Using 
right-handed electrons one could study a pure triple-gauge-coupling process, but this would re- 
quire longitudinally polarized electron beams, the prospect of which looks rather unfavourable 
for LEP2. Furthermore, for all energies the cross-section for right-handed electrons is sup- 
pressed by two orders of magnitude compared with the dominant left-handed mode, mainly 
because there is no t-channel contribution. Therefore essentially only the latter can be investi- 
gated at LEP2. With transverse beam polarization, however, one could obtain information on 
the right-handed matrix element via its interference with the left-handed one. 

As stated in section |2T^, for on-shell W-pair production and unpolarized electrons and 
positrons there are nine independent helicity-matrix elements or six if CP is conserved. These 
yield nine or six independent observables. Taking into account the decay of the W bosons there 
are many more observables: 81 or, if CP is conserved, 36 products of the various helicity-matrix 

'^This holds for arbitrary CP-conserving s-channel contributions in the limit of vanishing electron mass 
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[GeV] 


C"Born 


^Born 


^Born 


161.0 


3.812 


7.622 


0.002 


175.0 


15.959 


31.716 


0.202 


184.0 


17.427 


34.567 


0.287 


190.0 


17.762 


35.203 


0.321 


205.0 


17.609 


34.867 


0.350 



Table 2: Integrated lowest-order cross-section in pb for different polarizations of the electrons 
and different centre-of-mass energies. 



elements [ITJ]. Because the V—A structure of the W decays is well established, these observables 
can be extracted in a model-independent way from the data for the five-fold differential cross- 
section d(j/(dcos 6* dcos 6*1 d0i dcos6'2 d02), where 9i^2 and 0i^2 represent the decay angles of the 



two W bosons. These observables may serve to put limits on anomalous couplings |[Ty, [T^| 



Note that the five-fold differential cross-section requires the identification of the charge of at 
least one of the decaying W bosons, as otherwise the information on the sign of cos 6 would 
be lost. This is possible for hadronic-leptonic or leptonic-leptonic events. If no charge iden- 
tification is possible only production-forward-backward-symmetric observables would be left. 
As argued above, in those observables for instance the s-channel contribution is suppressed by 
P"^ in the threshold region^. Consequently, if one cannot identify the jet charges the purely 
hadronic events will not be of much use for studies of the gauge couplings at LEP2. Of course 
this does not concern the W-mass determination as it does not rely on s-channel contributions. 



2.3 Radiative corrections 

As has been argued in section the SM theoretical predictions for W-pair production should 
have an uncertainty of about 0.5% (2% at ^/s = 161 GeV) in order to obtain reasonable 
limits on the structure of the gauge-boson self-couplings, and to determine the W-boson mass 
with the envisaged precision of roughly 40-50 MeV. In this context radiative corrections are 
indispensable. 

Much effort has been made in recent years to obtain such precise theoretical predictions 
for W-pair production. In the following we discuss the existing results for the virtual and real 
electroweak corrections in the on-shell case. In addition we discuss the quality of an improved 
Born approximation (IB A) that contains all familiar, LEPl-like leading corrections for the W- 
pair production cross-section at LEP2 energies. Such a discussion is in particular relevant as 
the present off-shell LEP2 Monte Carlos often make use of such an approximation. 

•''if CP is violated there exists an anomalous gauge-boson coupling that does not yield a suppressed s-channel 
contribution |||, |ll| . 
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2.3.1 Radiative electroweak 0{a) corrections 



The 0{a) radiative corrections can be naturally divided into three classes, the virtual, soft- 
photonic, and hard-photonic contributions. Since the process e"'"e~ W"'"W~ involves the 
charged current in lowest order, the corresponding radiative corrections cannot be separated 
on the basis of Feynman diagrams into electromagnetic and weak contributions in a gauge- 
invariant way. Like we have already observed in equation ([T5|) , SU{2) gauge invariance requires 
an interplay between the 7- and Z-exchange diagrams. 

The complete radiative electroweak 0{a) corrections have been calculated independently 
by two groups [jl4|, |15|0. For the unpolarized case they have been checked and found to agree 
within a couple of per- mil (i.e. within the integration error). 

The virtual corrections contain infra-red (IR) divergences, which result from virtual photons 
exchanged between external charged particles. They are compensated for by adding the cross- 
section for the process e"'"e~ — > W^W~7. If the energy of the emitted photon is small compared 
with the detector resolution (soft photons), this process cannot be distinguished experimentally 
from the non-radiative W-pair production process. In practical experiments the soft-photon 
approximation is in general not sufficient and one has to include the radiation of hard photons, 
too. When adding these real-photon effects to the contribution of the virtual corrections, not 
only the IR singularities but also the large Sudakov double logarithms \og^{s/ml) drop out. 

Still there are various sources of potentially large 0{a) corrections left at LEP2 energies. 
First of all there are large QED corrections of the form {a/7T)log{Q'^/ml) with 3> ml, 
originating from collinear photon radiation off the electron or positron (see appendix ^). They 
form a gauge-invariant subset of QED corrections and amount to roughly 6% at LEP2, not 
taking into account possible enhancements from the corresponding coefficients (like, e.g., close 
to thresholds). 

From the renormalization two sets of potentially large fermionic (formally weak) corrections 
arise. The first set is associated with the charge renormalization at zero momentum transfer, 
where the relevant scale is set by the fermion masses entering the vacuum polarization. In high- 
energy experiments, however, the running charge should be evaluated at scales much larger than 
the masses of the light fermions fe = {e, /i, r, u, d, c, s, b}. This leads to large logarithmic ('mass 
singular') contributions of the form [a / tt) loglQ"^ / m'j^) with ^ m^^, which can amount to 
a shift in a of 8% at LEP2 energies. Related to the top quark, corrections oc /M^ will 
occur. They show up as universal corrections via the renormalization of the W and Z masses 
(or equivalently s^) if the corresponding renormalization scales are small compared with the 
mass splitting in the (t,b) isospin doublet. 

Finally the long-range electromagnetic interaction between the slowly moving W bosons 



leads to the so-called Coulomb singularity [Q. This singularity yields an 0{a) correction factor 



air/ (2/?), resulting in an 0(a)-corrected cross-section that does not vanish at threshold for left- 

^The process e+e^ W+W^7 has also been calculated in j2^. 
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Mh [GeV] = 


300 


300 


300 


60 


300 


1000 


«s = 


0.117 


0.123 


0.129 


0.123 


0.123 


0.123 


mt [GeV] = 


164.80 


165.26 


165.73 


148.14 


165.26 


181.20 


[GeV] 


^[ 


pb] 


161.0 


2.473 


2.472 


2.470 


2.514 


2.472 


2.460 


175.0 


14.471 


14.465 


14.459 


14.548 


14.465 


14.422 


184.0 


16.619 


16.613 


16.607 


16.671 


16.613 


16.568 


190.0 


17.263 


17.257 


17.251 


17.301 


17.257 


17.213 


205.0 


17.683 


17.677 


17.671 


17.698 


17.677 


17.637 



Table 3: Integrated unpolarized cross-section including radiative electroweak 0{a) corrections 
for different values for the Higgs-boson mass and as, at various centre-of-mass energies. The 
theoretical error in mt is roughly 1-2 GeV. 

handed electrons. The right-handed cross-section remains suppressed by at least This 
correction factor exhibits the fact that the free-particle approximation is inadequate near the 
W-pair production threshold in the presence of the long-range Coulomb interaction. We want 
to stress here that the Coulomb singularity is changed substantially by effects that effectively 
truncate the range of the interaction, like the off-shellness and the decay of the W bosons. This 



will be treated in section 3.4.2 and will lead to the conclusion that higher-order Coulombic 



corrections to the total cross-section are not important. 

The sensitivity of the total unpolarized cross-section to and the unknown mass of the 
Higgs boson is illustrated in table |^. The dependence on a^, originating from the calculation 
of mt, is completely negligible ( < 0.1%). Compared with the lowest-order cross-sections of 
table 0, a variation of Mh between 60 and 1000 GeV, however, influences the total cross-section 
by around 0.5% at the three highest energy points, by about 0.8% at ^/s = 175 GeV, and by 
1.4% at threshold. These numbers are lowered by about 0.5% if the radiative corrections are 
calculated in the so-called representation, which absorbs the universal mt and Mh effects 
present in the W wave- function factors (see section p. 4. 2] ). In this representation a variation 
of Mh between 300 and 1000 GeV has a negligible impact on the cross-section. Keeping in 
mind that we would like to reach a theoretical accuracy of 2% (0.5%) at ^/s = 161 GeV (i/i > 
175 GeV), it should be clear that the Higgs- mass dependence constitutes a major uncertainty. 



A more detailed investigation [^ revealed that this Higgs-boson-mass dependence is the result 
of the Yukawa interaction between the two slowly moving W bosons (mediated by the Higgs 
boson). As such the effect is largest close to threshold and for the lightest Higgs masses, which 
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a 
[pb] 




160 170 180 190 200 

V~s [GeV] 

Figure 2: Unpolarized total cross-section in the threshold region. The dotted curve corresponds 
to Born, the dashed one to G^-Born, and the solid one to the results including full radiative 
electroweak 0{a) corrections. 

yield the shortest range of the interaction. For instance, upon increasing the lower Higgs-mass 
bound from 60 to 90 GeV, the resulting uncertainty is reduced by roughly a factor of twoQ. 

In figure |^ we display the influence of the full 0{q) corrections on the total unpolarized 
cross-section near the W-pair production threshold. 

As has been stated before, it should be noted that, unlike at LEPl, the presence of the 
charged current at lowest order in e"'"e~ W^W~ complicates the separation between QED 
and weak corrections. So, in order to have a cleaner look at the QED corrections we present 
in addition to the Born cross-section also the lowest-order cross-section in the so-called 
representation, i.e. the Born cross-section with a replaced by \/2G^M^s\^/'k^ which already 
contains an important part of the leading weak effects discussed above]] (see also section p.3.3| ). 
We will refer to this cross-section as G^-Born. As a result of the steep drop with decreasing 
centre-of-mass energy of the W-pair cross-section close to threshold, large and predominantly 

''At the start of the 161 GeV run no significant change in the Mh bound is expected. So, only after the 
higher-energy LEP2 runs have taken place the improved knowledge on A/h can be used for an a posteriori 
reduction of the Mh dependence of the 161 GeV run. 

*As the t channel is dominant at LEP2 energies, the G^-Born describes the leading weak corrections reason- 
ably well for the default Higgs mass Mn^SOOGeV. 
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Figure 3: Unpolarized differential cross-section da/{dcos9^) at y/s = 190 GeV. The same 
signature as adopted in figure ||. 



soft QED effects can be observed. Compared with G^-Born the 0{a) corrections amount to 
more than —25% in the direct vicinity of the threshold (i/i < 165 GeV), apart from the region 
very close to threshold where the positive Coulomb- singularity contribution takes over (see 
close-up in figure D, and still -17% (-11%) at = 175 GeV (190 GeV). The size of these 
effects necessitates the inclusion of higher-order QED corrections in order to end up with an 
acceptable theoretical uncertainty. A discussion of these higher-order QED corrections can be 



found in section |2.3.3| . The finite W width has a drastic impact on the effects related to the 
Coulomb singularity (see section ^.4.2] ), but the large soft QED effects are merely smoothened 
and stay sizeable. 

In figure |^ the effect of the 0{a) corrections on the unpolarized differential cross-section 
dcr/ (dcos0+) is shown for ^/s = 190 GeV. Here ^+ stands for the polar angle of the W"*" boson 
with respect to the incoming positron^. Apart from the expected normalization effects that 
are already observed in the total cross- sect ion, also a distortion of the distribution occurs. As 
this is exactly the type of signature one might expect from anomalous gauge-boson couplings, 
this underlines again the importance of having a profound knowledge of the SM radiative 
corrections. The origin of the distortion can be traced back to hard initial-state photonic 
corrections. Hard-photon emissions boost the centre-of-mass system of the W"^W~ pair. As a 

■'In the presence of hard-photon radiation in general 9+ 
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result of that, events that are forward in the centre-of-mass system of the produced W bosons 
can show up as backward events in the LAB system and vice versa. Since the cross-section in 
the backward direction is substantially lower than in the forward direction, the net effect of 
this redistribution (migration) of events will be a distortion of the differential distribution with 
respect to the lowest-order one. Of course these boost effects are much more pronounced at 
high energies 0. 



2.3.2 Approximations in the LEP2 region 

The complete analytic results for the electroweak 0{a) corrections are very lengthy and com- 
plicated, resulting in huge and rather slow computer programs. Moreover, the formulae are 
completely untransparent. In view of this, simple approximative expressions are desirable. 
Apart from providing fast computer programs, which are useful for many applications, simple 
transparent formulae should reveal the physical content and the origin of the dominant radia- 
tive corrections. Furthermore, if these approximative expressions represent the exact corrections 
adequately, one might consider implementing them in the existing LEP2 Monte Carlos. 

Owing to the lack of a calculation of the complete 0{a) corrections, the present LEP2 Monte 
Carlos for off-shell W-pair production include only the known leading universal corrections. In 
order to assess the theoretical uncertainty inherent in this approach, the on-shell case can be 
used as guideline. The size of the non- leading 0{a) corrections in this case should provide a 
reasonable estimate for the corresponding left-out non-leading corrections in the off-shell case. 

We start out by investigating the structure of the matrix element for e^e~ W"'"W~. 
Whereas it involves only three different tensor structures in lowest order, at 0{a) twelve in- 
dependent tensor structures occur, each of which is associated with an independent invariant 
function, which can be considered as an s- and t-dependent effective coupling. The dominant 
radiative corrections, as e.g. those that are related to UV, IR, or mass singularities, in gen- 
eral have factorization properties and are at 0{a) restricted to those invariant functions that 
appear already at lowest order. Therefore the contributions of the other invariant functions 
should be relatively small. Indeed, a numerical analysis reveals that in a suitably chosen rep- 
resentation for the basic set of independent matrix elements only the three Born-like invariant 
functions plus one extra right-handed piece, related to Aif, are relevant for a sufficiently good 



approximation This results in the following approximation for the matrix element: 

m:^^ = m^jFp + m:^f^, (19) 



with A4j and A4q defined in equation ([T5|). The term involving Fj is needed for right-handed 
electrons because of contributions cx (3"^ cos 6, originating from the interference of A^Bom with 
particular 1-loop matrix elements, which are not present at lowest order. Neglecting the other 



invariant functions in the basis given in ||32| introduces errors well below the per-cent level 
(see tabled). This is of course only true for observables where A^app is not suppressed or 
absent. All this demonstrates that improved Born approximations are possible. Note, however. 



17 



that in contrast to the situation at LEPl the invariant functions F/g are both energy- and 
angular-dependent . 

In order to construct an improved Born approximation (IBA) one has to specify simple ex- 
pressions for the invariant functions Fpq, which reproduce the corresponding exact expressions 
with sufficient accuracy. In the LEP2 energy region the following expressions can be used as a 
reasonable ansatz [ P^ 



Airais) + Ana '^(1-0 
4/? 



2\2 



(20) 



The terms containing and a{s) incorporate all the leading universal corrections associ- 
ated with the running of a and the corrections oc am^/M^ associated with the p parameter 
(see also section p.3.3| ). As these are linked to the renormalization of the electric charge at zero 
momentum transfer and of the weak mixing angle, they only contribute to the structures al- 
ready present at lowest order. The corresponding leading 0{a) contributions can be recovered 
from equation (BO) by substituting 



a 



» ^ a[l + Aa(s)], 2v^G^M 



with 



w 



m 



Aira 



2 / ' 



'W 



Ap 



3am^ 



(21) 



(22) 



The 1//3 term describes the effect of the Coulomb singularity, which for f3 <^ 1 yields a 
simple correction factor to the lowest-order cross-section: 



na 

OCTCoul = 0"Born- 



(23) 



The factor (1 — is introduced by hand to switch off the Coulomb singularity fast enough. 
In this way one avoids corrections that are too large away from threshold, where the Coulomb 
singularity should not play a role anymore. As has been studied in heavy-mass contribu- 
tions oc log(mt) and log(MH) that are not covered by equation (pO]) have a negligible impact 
on the approximation. Since the above IBA analysis has been performed for the default Higgs- 
boson mass Mh = 300 GeV, the large light- Higgs-boson corrections close to threshold, shown 
in table |^, are absent. By adding a simple approximation for these corrections to the IBA, 
however, the full Mh dependence of the exact 0{a) corrections can be reproduced. 

In addition to the contributions described so far, one has to include the leading logarithmic 
QED corrections. These can be calculated using the structure-function method as described in 
appendix^. They comprise all contributions oc (a/vr) \og{m1/Q'^). In the following numerical 
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analysis the scale = s has been used and, in order to extract the effect of the non-leading 
virtual corrections, hard-photon radiation is left out. 

As we want to compare the approximation with an 0{a) calculation we do not use the 
square of the matrix element defined in equations ( jlQ]) and ( pO]) for the numerical analysis, 
but only the square of the terms involving and a{s) and the interference of these terms 
with the others. Moreover, in this interference and a{s) are replaced by the corresponding 
lowest-order expressions in terms of and a(0). Nevertheless, the so-obtained approximation 
still contains higher-order contributions through the squares of a{s) and G^. To allow for a 
meaningful comparison these have been included in the numbers for the exact one-loop results 
given in table § in the same way (see |32] for more details). 



In table § we show the difference between the exact and approximated virtual and soft 0{q) 
radiative corrections to the total and differential cross-section for right-handed, left-handed, and 
unpolarized electrons. The positrons are assumed to be unpolarized. In the LEP2 energy region 
the relative difference between the exact result and the approximation can be as large as 1-2% 
for left-handed or unpolarized electrons, and reach 3-8% for right-handed electrons. As there 
is no obvious reason why these remaining non-leading corrections should be smaller in the case 
of off-shell W bosons or e+e~ 4/, their inclusion in the LEP2 data analysis seems to be 
unavoidable. 



2.3.3 Higher-order corrections 

In section |2.3.1| we have encountered large 0{a) corrections. A short description of the way to 
include the corresponding higher-order corrections is hence in place. 

In two distinct ways the higher-order corrections enter the calculation of the distributions 
for W-pair production. First of all there is the calculation of mt from Mw and the input 
parameters, either to be used in the calculation of the W-pair RC's or to be confronted with 
the Tevatron data. As mt enters the relation between a, G^, Mz, and Mw, resulting from 
the muon decay width, at the 1-loop level, the highest precision possible is required for this 
relation. As this relation relies on calculations performed at ~ for the muon decay width 
and at the subtraction points = M^, M| for the renormalization procedure, we can use the 
state-of-the-art calculation developed for the LEPl analysis [|]. This yields 

"^^v^'^- '''' 

Usually Pc is written in the form 
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unpolarized 


right-handed 


left-handed 


v/i = 161 GeV 


total 


1.45 


0.00 




-1.56 


-0.01 


1.45 


0.00 


10° 


1.63 


0.00 




4.41 


0.00 


1.63 


0.00 


90° 


1.44 


0.00 




-1.57 


-0.01 


1.44 


0.00 


170° 


1.26 


0.00 




-7.52 


0.00 


1.26 


0.00 






V 




" = 165 GeV 






total 


1.27 


0.02 




-2.09 


-0.01 


1.28 


0.02 


10° 


1.67 


0.00 




0.49 


0.00 


1.67 


0.00 


90° 


1.17 


0.02 




-2.09 


-0.02 


1.18 


0.02 


170° 


0.75 


0.00 




-4.64 


0.00 


0.77 


0.00 






V 




"= 175 GeV 






total 


1.26 


0.03 




-2.58 


-0.02 


1.28 


0.03 


10° 


1.71 


0.00 




0.18 


0.00 


1.71 


0.00 


90° 


1.03 


0.05 




-2.59 


-0.03 


1.06 


0.05 


170° 


0.59 


0.00 




-5.30 


0.00 


0.69 


0.00 






V 




"= 184 GeV 






total 


1.02 


0.04 




-2.80 


-0.03 


1.06 


0.04 


10° 


1.57 


0.00 




2.17 


-0.01 


1.57 


0.00 


90° 


0.67 


0.08 




-2.82 


-0.05 


0.72 


0.08 


170° 


0.10 


0.00 




-7.72 


0.01 


0.32 


0.00 


v/i = 190 GeV 


total 


1.24 


0.03 




-2.91 


-0.04 


1.28 


0.03 


10° 


1.67 


0.00 




0.59 


-0.01 


1.67 


0.00 


90° 


0.95 


0.06 




-2.92 


-0.06 


1.01 


0.06 


170° 


0.58 


0.00 




-6.32 


0.00 


0.83 


0.00 


v/i = 205 GeV 


total 


1.60 


0.00 




-3.11 


-0.09 


1.65 


0.00 


10° 


1.77 


0.00 




-1.68 


-0.01 


1.77 


0.00 


90° 


1.55 


0.00 




-3.14 


-0.12 


1.64 


0.00 


170° 


1.61 


0.00 




-4.37 


0.00 


1.94 


0.00 



Table 4: Quality of the IBA (first column per polarization) and the form-factor approximation 
(second column per polarization), given in per cent relative to Born. Results are given for 
the total cross-section (total) and the differential cross-section at 10, 90, and 170°. 
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where Ar contains all the one-loop corrections to the muon decay width with the inclusion 
and the proper arrangement of the higher-order terms. Next we subdivide pc as introduced in 



equation ( pSj) into a leading term, Ar^, and remainder terms, Ar^em, as follows: 



Ar 1 



Ar^ - Ar^em 



(26) 



Aa) ( 1 + =fAp 



Ar, 



'w 



with 



Ar, 



Ar" + Ar""= + ^Ap^ - Aa 



(27) 



This contains all the terms known at present: the complete one-loop 0{a) corrections Ar" 
(two-, three-, four-point functions) and complete two-loop 0{aas) insertions Ar""" to two- 
point functions, from which the leading 0{a) and 0{aas) terms are subtracted. The Born 
version of equation (^31) , i.e. pc = 1, relates directly to a, Mw, and Mz, hence leaving no 
room for using as additional independent input parameter. However, having introduced 
radiative corrections to the Born relation we can solve equation (|^) in terms of the top quark 
mass and, in turn, this particular value for rrit will be used throughout the rest of the calculation. 
One should notice that this procedure is not free of ambiguities since Ar is also a function of 
Mh, but with much smaller sensitivity due to the well known screening. So in the end we 
should also allow for some variation in Mh. Moreover, Ar is also a function of through 
higher-order corrections, for instance 0{aas), and, as a consequence, also some variation in the 
strong coupling constant should be allowed. 



The leading terms Aa and Apx appearing in equation (p7|) are given by 



Apx = Ap" + Ap""» + X 



3a 



Aa 



ais] 



167rs^Cw 
a 



2 

Ml 



1-- 1 + — 
3 I 3 



7r2\ a,(m2) 



TT 



a 



{2i 



where T}^^{s) is the renormalized 0{a) transverse photon self-energy originating from fermion 
loops, excluding top-quark loops. 

The term X in equation ( ^8]) is a next-to-leading order term, whose proper treatment is 
rather important: 



X = Re 



nz(M|) 

M| 



nw(M^) 



Hoop 



MS 



Ap^ 



(29) 



where IIv denotes the unrenormalized transverse self-energy of the V gauge boson. The ultra- 
violet (UV) divergences are removed according to the MS renormalization scheme with p = Mz- 
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In contrast to Apx, the leading contribution Apx appearing in equation (p6|) is normalized 



by rather than by a/(s^M-^), as is required by the resummation proposed in |^ 



Apx = Ap" + Ap" + Ap^ 



N'cxt 



1 +xt Ap 



(2) / "^t 



+ Ap° 

2 \ 



+ X 



Ml, 



+ Ci 



TT 



+ C2 



where Nq 



3 and 



X 



Gf, ml 



(30) 
(31) 

- (32) 

The coefficients ci and C2 describe the ffist- and second-order QCD corrections for the leading 
Xt contribution to Ap, calculated in |^ Correspondingly: 



X. 




Cl = 

C2 = -vr^ (2.564571 - 0.180981 Uf) , 

with n/ the total number of flavours {nj = 6). The function Ap'^^)(m^/M^' 
leading two-loop electroweak xt correction to Ap, calculated ffist in the Mh = 
in ^% and later in |3^ for an arbitrary relation between Mh and mt- 



(33) 
(34) 

describes the 
approximation 



It should be noted that the higher-order radiative corrections discussed above are usually 
calculated in the limit of a heavy top mass, i.e. usually only the leading part of the corrections 
is under control. This often raises the question of what effect can be estimated from the sub- 
leading terms, since m^ ~ 2Mz. As far as QCD 0{as) and ^(a^) corrections are concerned 
the sub-leading terms are well under control. In the pure electroweak sector, however, there 
has been, so far, no calculation of higher-order sub-leading terms at an arbitrary scale. The 
only available piece of calculation concerns the p parameter at = B^], therefore relevant 



only for t'^-e scattering. If however one is willing to extrapolate the Q = result to a higher 
scale, by assuming that the ratio of leading to next-to-leading corrections is representative for 
the corresponding theoretical uncertainty, then the answer is next-to-leading ^ leading. 

Even before considering the actual impact of the higher-order terms on mt, we should men- 
tion at this point that the way in which the non- leading terms can be treated and the exact 
form of the leading-remainder splitting give rise to several possible options in the actual im- 
plementation of radiative corrections. This in turn becomes a source of theoretical uncertainty. 
For instance, for Ar we can introduce the decomposition into leading and remainder. Since we 
know how to proceed with all objects in the leading approximation, the only ambiguity is due 
to the treatment of the remainders. Clearly, after the splitting Ar = Ar^ -|- Arrem there are in 
principle several possible ways of handling the remainder: 

1 Arrcm \ 1 + Arrem 1 



Ar 



Ar^ — Arr 



Ar, 



1 + 



1 - Ar, 



1 - Ar, 



Ar, 



+ Ar„ 



(35) 
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Actually, these options differ among themselves, but the difference can be related to the choice 
of the scale in the remainder term. A complete evaluation of the sub-leading (9(G^M|m^) 
corrections would greatly reduce the associated uncertainty. In conclusion, we observe that a 
natural and familiar language for the basic ingredients of the physical observables is that of 
effective couplings. However, it should be stressed that the above formulae belong to a specific 
realization of this language and other realizations could also be used. 

Now we can assess the influence of the higher-order corrections on the calculation of rrit 
from equation (0). Using for instance Mw = 80.26 GeV and Mz = 91.1884 GeV we find 

mt = 165li^(MH,«.)GeV, (36) 

the central value of which corresponds to Mh = 300 GeV and as(M^) = 0.123. The errors are 
derived by varying Mh and in the range 60 GeV < Mh < 1 TeV and = 0.123 ± 0.006. It 
should be noted in this respect that the total variation induced by Mh alone (at = 0.123) is 
about 33 GeV. The following is observed for the higher-order corrections: 



• by neglecting the (9(a^) term in Ar we find for the same input parameters (and Mh = 
300 GeV, as = 0.123) a shift in rrit of -1.9 GeV. 

• If instead we switch off the 0(0:0;^) correction the corresponding shift will be —1.5 GeV. 

• If finally both the 0{aas) and 0{aa1) corrections are neglected we find a shift in rrit of 
— 10.5 GeV. Here by 0{aas) the full result is implied and not only the leading part of it. 



Based on the above observations, the remaining theoretical uncertainty in the calculation of 
rrit from missing higher-order corrections and sub-leading 0(q;^) corrections to Ar is estimated 
to be roughly 1-2 GeV. 

The second way the higher-order corrections enter the calculations for W-pair production is 
through the process itself, so through self-energies, vertices, etc. The known weak higher-order 
effects comprise the running of a [see equation (^)] and the complete 0{aas) corrections to 
the gauge-boson self-energies |35|. Other higher-order calculations, as those for Ap at 0{a^) or 
0{aal), have been performed in the limit where rut is the largest scale. The QCD corrections 
associated with the gauge-boson vertex corrections have, to our knowledge, not been calculated 
yet, but they are at most logarithmic in the top mass. If one takes into account the leading 
weak effects by replacing e^/(2s^) by 2^/2G^M^ in the A^/ part of equation ([T5| ) and a by 
a{s) in the Aiq part, the remaining unknown higher-order weak effects are expected to be 
negligible compared with the required theoretical accuracy. 

As pointed out in the previous subsection, the virtual and real corrections reveal the presence 
of large logarithmic QED effects of the form aL/vr = (a/vr) log{Q'^/ml) with ^ ml. They 
arise when photons or light fermions are radiated off in the direction of incoming or outgoing 
light particles. In the case of W-pair production the only terms of this sort originate from initial- 
state photon emission. Radiation of photons from the final-state W bosons can only lead to 
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JDUi ii 


+ h.o.t. for Q'^ = s 


+ h.o.t. for = s 


frpVi 

[vjc V J 


-J- n(r\\ 


0{a'^) LL 


exp. LL 


0\a^) LL 


exp. LL 


161.0 


2.472 ± 0.001 


3.103 


3.003 


3.095 


2.998 


165.0 


8.581 ± 0.003 


9.079 


9.049 


9.061 


9.033 


170.0 


12.270 ± 0.004 


12.583 


12.585 


12.567 


12.568 


175.0 


14.465 ± 0.004 


14.654 


14.670 


14.642 


14.656 


184.0 


16.613 ± 0.005 


16.668 


16.693 


16.663 


16.686 


190.0 


17.257 ± 0.006 


17.259 


17.286 


17.259 


17.283 


205.0 


17.677 ± 0.006 


17.613 


17.638 


17.620 


17.641 



Table 5: Unpolarized total cross-section, given in pb, including radiative electroweak 0[a) 
corrections and higher-order terms in the leading-log approximation. These higher-order terms 
are given with and without soft-photon exponentiation for two different 'natural' scale choices. 

sizeable corrections if the energy of the W bosons is much larger than their mass. The leading 
large logarithmic corrections can be calculated by using the so-called structure-function method 
]16| in leading-log (LL) approximation, i.e. only taking along the terms oc (aL/vr)". This 



procedure also allows the inclusion of soft-photon effects to all orders by means of exponentiation 
and is discussed in detail in appendix ^. 

In table ^ the higher-order effects related to the large logarithmic QED corrections are 
displayed. The 0(q;^) LL entry contains, in addition to the full 0{a) result, the contribution 
from 0[o?^ LL corrections using o"o = crBorn in the convolution (]68|). The exp. LL entry 
contains on top of that the exponentiation of soft-photon effects. Compared with the full 
0{a) results we observe large 0(a;^) LL effects near threshold, e.g. 25% at ^/s = 161 GeV, 
and moderate ones when sufficiently far above threshold, i.e. < 1% for energies above roughly 
175 GeV. The additional soft-photon exponentiation is only sizeable near threshold, e.g. —4% 
at ^/s = 161 GeV. Note that finite- W- width effects will smoothen the threshold behaviour 
and hence reduce the size of the 0{a'^) LL effects, nevertheless they will stay sizeable near 
threshold. The dependence of the higher-order LL corrections [beyond on the scale 

choice is negligible, since all natural scales are roughly equal close to threshold. When 
comparing the popular scale choice Q"^ = s with = s(l — + P), motivated by the 

behaviour of the total cross-section near threshold and at high energies 0, the differences are 
at the 0.1% level (0.2-0.3% at ^ = 161 GeV). 
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An additional improvement of the theoretical predictions can be obtained by using an im- 
proved Born cross-section in the convolution (|68D , taking into account corrections related to 
G^, Mh etcQ 



2.4 The width of the W boson 

Evidently the width of the W boson is a crucial ingredient for the (off-shell) W-pair production 
cross-section, especially in the threshold region. Moreover, the branching ratios enter the cross- 
sections for definite fermions in the final state. As at present the width of the W boson is 
experimentally poorly known, we need a precise theoretical calculation instead in order to 
obtain adequate theoretical predictions for off-shell W-pair production. 



2.4.1 The W-boson width in lowest order 

The W width is dominated by decays into fermion-antifermion pairs. In lowest order the partial 
width for the decay of a W boson into two fermions with masses mj. and mf denote the 
generation index and /, /' stand for u, d or u, 1) is given by 



_,Born 



6 2^2^ ' 



m 



2M^ 



(37) 



For leptonic decays the mixing matrix is diagonal (Vij = 6ij) and the colour factor Nq equals 
one. For decays into quarks there is a non-trivial quark mixing matrix and = 3. For the 



quark mixing matrix we have used Si2 = 0.221, S23 = 0.04, and S13 = 0.004 [^. The total 
width is obtained as a sum over the partial fermionic decay widths with rrif - + mf < Mw 



_^Born 
W 



E Born ^ poui" 



Born 



(3^ 



As the quark masses are small compared with Mw, the fermion-mass effects are small for the 
W decay. If we neglect them we obtain the simple formulae 



Born _ / a Mw 
i \\r f. f — -/ V( 



W/i/' 



6 2s^ 



Born 
I cy 



3a M 



w 



2 2sw 



(39) 



'^As the finite decay width of the W bosons wiil have a substantial impact on the Coulomb singularity, a LL 
analysis involving this Coulomb singularity only makes sense for off-shell W bosons. 
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Mh [GeV] 


300 


300 


300 


60 


300 


1000 


as 


0.117 


0.123 


0.129 


0.123 


0.123 


0.123 


mt [GeV] 


164.80 


165.26 


165.73 


148.14 


165.26 


181.20 


r7° [GeV] 


1.9490 


1.9490 


1.9490 


1.9490 


1.9490 


1.9490 


r7° [GeV] 


2.0354 


2.0354 


2.0354 


2.0354 


2.0354 


2.0354 


Tw [GeV] 


2.0642 


2.0663 


2.0684 


2.0681 


2.0663 


2.0639 


Tw [GeV] 


2.0791 


2.0817 


2.0844 


2.0813 


2.0817 


2.0819 




0.03416 


0.03398 


0.03380 


0.03491 


0.03398 


0.03275 




-0.00347 


-0.00347 


-0.00347 


-0.00369 


-0.00347 


-0.00341 


^QCD 


0.02495 


0.02623 


0.02751 


0.02623 


0.02623 


0.02623 



Table 6: Higgs-mass and dependence of the W-boson width 



2.4.2 Higher-order corrections to the W-boson width 



The electroweak and QCD radiative corrections for decays into massless fermions (m/ -C Mw) 
have been calculated in [|l3]-||l9l- The full one-loop electroweak and QCD corrections, together 



with the complete photonic and gluonic bremsstrahlung, were evaluated for arbitrary fermion 
pO|] . The various calculations are in good agreement. The relative corrections to the 



masses m 



total W-boson decay width are given in table ^ The electroweak corrections in the on-shell 



scheme S™ = Tt^/T^^^ — 1 — S^^^ are predominantly originating from the running of a and 
the corrections oc am\/M^ to the p parameter. These corrections can be easily accounted for 
by parametrizing the lowest-order width in terms of and instead of a and s\f. The 
width in this parametrization Fw is related to the width in the on-shell parametrization 
Tw by 



w — J- w 



w — ^ X ~ ■'^w 



1 - Ar 1 - Ar 

= T^'\l + + - Ari-'°°P) ^ r7"(l + + (40) 

As can be seen from table |] the electroweak corrections with respect to the parametrization 
with G^, 5''™, depend in a negligible way on Mh, and remain below 0.5% for the total width. 
The QCD corrections 5'^'^^ are practically constant and equal to 2as(M^)/(37r), their value for 
zero fermion masses. For the numerical evaluation we use as(M-^) = 0.123 (i.e. equal to the 
default input value). The difference between Fw and the more precise F-w is caused by missing 
higher-order terms related to Ar. 



We obtain the following improved Born approximation for the total and partial widths [^ 
Fw,.,L- ^ F 



TBorn 



TT 
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Born 
nif 7^ 


complete 
ruf 7^ 


complete 
rUf = 


IBA 

(m/ = 0) 


Branching 
ratio 


r(W ^ ez/e) 


0.2262 


0.2255 


0.2255 


0.2262 


0.1083 


r(w -> /iz/^) 


0.2262 


0.2255 


0.2255 


0.2262 


0.1083 


r(W ^ TUr) 


0.2261 


0.2253 


0.2255 


0.2262 


0.1082 


r(W lep.) 


0.6785 


0.6763 


0.6765 


0.6787 


0.3249 


r(w ^ ud) 


0.6455 


0.6684 


0.6684 


0.6708 


0.3211 


r(w -> us) X 10 


0.3315 


0.3432 


0.3432 


0.3444 


0.0165 


r(W ^ ub) X 10^ 


0.1080 


0.1122 


0.1124 


0.1128 


0.000005 


r(W cd) X 10 


0.3312 


0.3431 


0.3432 


0.3444 


0.0165 


r(w ^ cs) 


0.6441 


0.6672 


0.6673 


0.6697 


0.3205 


r(W cb) X 10^ 


0.1080 


0.1121 


0.1124 


0.1128 


0.0005 


r(W had.) 


1.3569 


1.4054 


1.4055 


1.4104 


0.6751 


r(W ^ all) 


2.0354 


2.0817 


2.0820 


2.0891 





Table 7: Partial and total W-decay widths Fw in different approximations given in GeV. 
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In table ^ we compare the improved Born approximation (IBA) for the partial and total 
widths, given by equation (|4l|), with the lowest-order widths, the widths including the complete 
first-order and leading higher-order corrections for finite fermion masses, and the same for 
vanishing fermion masses, all in the parametrization. The effects of the fermion masses, 
which are of the order m'j/M^, are below 0.3%. Consequently the exact numerical values 
for the masses of the external fermions are irrelevant. The IBA reproduces the exact results 
within 0.4% (0.6% for the decays into a b-quark). The branching ratios for the individual 
decay channels derived from equation (^I]), which depend only on and Vij, agree numerically 
within 0.1% with those obtained from the full one-loop results. 



3 OfF-Shell W-Pair Production 
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Figure 4: Lowest-order diagrams for e^e —>■ W+W 4/ 





du 


sc 








du 


43 


11 


20 


10 


10 


ei/e 


20 


20 


56 


18 


18 




10 


10 


18 




9 



Table 8: Number of Feynman diagrams for W-pair produced four-fermion final states. 
3.1 Lowest order: an introduction 



So far we have only considered the production of stable W bosons. This is, however, only an 
approximation and in particular in the threshold region it is not sufficient. Rather, one has to 
describe the W bosons as resonances, with a finite width so as to avoid singularities inside the 
physical phase space, and analyse their presence through their decay products: 



e+ + e- ^ W+ + W- ^/i + /2 + /3 + /4. 



(42) 



Process ( ^2]) involves two resonant W bosons (doubly-resonant) and can be viewed as a very 
natural first step beyond the on-shell limit. In lowest order this process is represented by the 
three Feynman diagrams shown in figure ^ However, the full four-fermion process does not only 
proceed through the three doubly-resonant diagrams. There are also contributions from other 
diagrams with the same initial and final states, but different intermediate states. Classifications 
of four-fermion production processes and of the contributing diagrams are given in [^, ^ . In 
table H we give the number of diagrams contributing for final states that can be reached by 
W-pair intermediate states. These so-called charged- current processes are sometimes referred 
to as CCn, with n denoting the number of contributing diagrams [e.g. CCS denotes process 
(|i2|)] . The simplest case (boldface numbers in table ^ is fully covered by doubly- and singly- 
resonant diagrams as given in figures | and |^. Additional graphs of the types shown in figure |] 
must be taken into account if electrons or electron-neutrinos are produced (roman numbers in 
table H). If the produced final state consists of particle-antiparticle pairs, the final state can 
also be obtained through intermediate Z-pair production, leading to extra Feynman diagrams 
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Figure 5: Example of a singly-resonant diagram. 



Figure 6: Examples of additional diagrams for final states with electrons and electron- neutrinos. 

{italic numbers in table |]). Finally, it should be mentioned that QCD diagrams, involving 
an intermediate gluon, contribute in the case of final states consisting of two quark-antiquark 
pairs. 



3.2 Semi-analytical approach 

We will now, in a first step of sophistication beyond on-shell W-pair production, introduce a 
finite W width and perform the CCS calculation. In the following a semi-analytical method will 



be emphasized, whereas Monte Carlo methods are treated elsewhere |^ . The starting point is 

= yds+pw(s+) J ds_pw(s-) fT^'^'(s;s+,s_), (43) 



with 5+ = k"^ and s_ = the virtualities of the internal W bosons. The Breit-Wigner densities 

vr \s± — + zMwrwl 

contain the finite width of the W boson, the coupling constants of its decay to fermions, and 
the corresponding branching ratio. In the limit of stable W bosons, the on-shell cross-section 
is recovered via 

Pw(s±) "^^^ - M^) X BR. (45) 
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The two-fold differential cross-section contains terms corresponding to the Feynman diagrams 
of figure | and their interferences. In the notation of it may be described by three terms, 



The couplings Cccs*'* contain the weak mixing angle, the Z-e vector and axial-vector coupling, 
the triple-gauge-boson couplings, and the s-channel Z and 7 propagators. The kinematical 
functions Gcc3^'^^ s-'^^ known analytically pT |: 



ns /j33 _|_ nst /j3/ , r,t nf f 
'-'CC3 ^CC3 ^ ^CC3 ^CC3 ^ '-'CCS ^CC3 



(46) 



1 

48 



A(s, s+, s_) + 12 s (s+ + s„) - 48 s+ s 
+ 24 s+ s_ (s — s+ — s_) £(s; s+, s_) 



^CC3 ('^j ■5+) '5- 



A(g,g+,g- ) p 
192 



A(s, s+, s_) + 12 (s s+ + s s_ + s+ s_ 



(47) 



(4J 



48 



— 24 s_ (s s_|_ + s s_ + s+ s_) C{s] s+, ■S-)! 



where 



A(s, s+, s_) = -|- + — 2 s s+ — 2 s s_ — 2 s+ s 
1 



£(s;s+,s_ 



A(s,s+,s_) 



log 



s_ + JA(s,s+,s. 



s_ - JA(s,s+,s_) 



(49) 

(50) 
(51) 



The numerical importance of the W-boson off-shellness in equation ( |45| ) is displayed in figure |^. 
Clearly, the effect of the finite W width is comparable to the one from universal initial-state 
radiation, defined in section |3.4.1| . 

The second level of sophistication of W-pair production is reached by performing calculations 
for complete sets of Feynman diagrams for specific final states. The simplest case, the CCll class 
of processes, corresponds to four-fermion final states without electrons or electron-neutrinos 
and without particle-antiparticle pairs. The corresponding gauge-invariant set of Feynman 
diagrams consists of the three doubly-resonant diagrams of figure ^ plus at most eight singly- 
resonant diagrams of the type shown in figure ^. This class of processes is still calculable in 
the semi-analytical approach and is of specific interest for the W-mass determination and the 
TGC studies. The corresponding cross-section is given by |?3| 



a 



cell. 



s) 



ds_i_ / ds 



A(s,s+,s_ 



15 

E 

fc=i 



ds+ ds_ 



(52) 
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with 

-. ^ — = CkGk{s,s+,S-) . (53) 

Initial- and final-state couplings, as well as intermediate W, Z, and 7 propagators are collected 
in the C^s. The kinematical details are contained in the QkS. Because of symmetry properties, 
under the exchanges of the three virtualities, only six of these kinematical functions are inde- 
pendent. Apart from the three kinematical functions of the CCS process, three more enter as a 



result of the singly-resonant diagrams. Explicit expressions can be found in |^ 



A similar semi-analytical analysis of the remaining four-fermion final states has not been 
performed so far, although this would be desirable. 
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3.3 Lowest order: the gauge- invar iance issue 



The discussion in the previous subsections has avoided the question of gauge invariance which 
arises when going from on-shell W-pair production to the off-shell case. There are two sources 
of gauge non-invariance one has to be aware of. 

First of all there is the issue of gauge non-invariance as a result of incomplete sets of 
contributions. Consider to this end again the process 



e+ + e- ->W+ + W- -./1 + /2 + /3 + /4, 



(54) 



represented by the three Feynman diagrams shown in figure ^ Forgetting about the finite W 
width the corresponding matrix element can be written in the following form: 

^2 



-MS 



Born 



2s^ 



w 



X 



X 



Ui7p uj- f 2 X 



X 



where 



2M| 



Ml 



v+ 



'•.+gP'' -YK + l'^^- 



I _ s 



(55) 



M 



(56) 



are directly related to the corresponding on-shell quantities J^j and A4q defined in equation ([T5|). 
Here k± are the momenta of the bosons as reconstructed from the decay products, and s 
and t the usual Mandelstam variables for e"'"e~ — W^W~. We have assumed that all fermion 
masses are negligible. In a renormalizable gauge such as a gauge one would otherwise have 
to consider in addition diagrams involving one or two unphysical charged Higgs bosons instead 
of the W bosons. 

Whereas for on-shell W bosons the contributions of Adj and Aiq are separately gauge- 
invariant, for off-shell W bosons not even A^Bom is gauge- invariant. This can be illustrated 
by considering the doubly-resonant diagrams in an axial gauge To this end we replace the 
conventional 't Hooft-Feynman gauge-fixing term for the W boson by —{n^W^){n'^Wj^). This 
modifies the W-boson propagators in the diagrams of figure ^ and, since this gauge-fixing term 
does not eliminate the mixing between the W-boson field and the corresponding unphysical 
Higgs field, leads to mixing propagators between those fields, giving rise to additional Feynman 
diagrams. Combining all relevant diagrams yields the following extra contribution to A^Born 
the axial gauge: 



Born 



X 



axial gauge 
MP 



Born 



{k't 



M^){n-kJ] 



't Hooft— Feynman gauge 

M^){n-k^ 



Mi 



9^2 



V2, 



■ Ui^p UJ^ V2 X 



e 



V2. 



■U3jaU_V4. (57) 
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Note that these gauge- dependent terms involve either a pole in (A;^ — M^) or (A;^ — M^) but 
not both, i.e. they are only singly- resonant. These terms are exactly cancelled by the gauge- 
dependent contributions of eight singly-resonant diagrams contributing to the same final state 
with the topology shown in figure |^. 

The sum of the doubly-resonant diagrams and the singly-resonant ones with the topology 
shown in figure ^ is in general gauge-invariant. This follows directly from the fact that for 
final states with four different fermions and no electrons or positrons those diagrams are the 
only ones that contribute. As a consequence the non-resonant diagrams are not needed to 
cancel the gauge- dependent terms of the doubly-resonant diagrams. Among the non-resonant 
diagrams and the singly-resonant ones that do not have the topology shown in figure |^ further 
gauge-invariant subsets can be identified by considering other four-fermion final states. 

Thus, in general all graphs that contribute to a given final state have to be taken into 
account and one is lead to consider the complete process e"'"e~ — > 4/, including all resonant 
and non-resonant graphs, in order to obtain a manifestly gauge-independent result. 

Simple estimates indicate that all non-doubly-resonant contributions are typically sup- 
pressed by a factor Fw/Mw ~ 2.5% on the cross-section level for each non-resonant W prop- 
agatoi||. This is confirmed by explicit calculations For instance, using covariant (i?^) 



gauges the universal non-doubly-resonant graphs that occur for all final states (see figure ||) 
contribute less than 0.15% to the total cross-section for 175 GeV < ^/s < 205 GeV, and 0.3% 
at 161 GeV. The contribution of the non-doubly-resonant t-channel photon-exchange graphs, 
which only occur when there are electrons in the final state, depends very much on the angular 
cut imposed on the outgoing electrons; for 10° they contribute at the per-cent level, e.g. ~ 4% 
at y/s = 190 GeV. Although all these non-doubly-resonant contributions are suppressed, they 
should nevertheless be taken into account to reach a theoretical accuracy of ~ 0.5%. Finally, 
the QCD graphs have been shown not to interfere with the electroweak graphs to any sizable 
extent, and can thus be computed using standard (Monte Carlo) QCD programs. 

Even when considering the complete set of graphs contributing to a given final state, there is 
still a more fundamental gauge-invariance problem to be solved. The resonant graphs discussed 
above involve poles at k"^ = These have to be cured by introducing the finite width 

in one way or another, while at the same time preserving gauge independence and, through 
a proper energy dependence, unitarity. In field theory, such widths arise naturally from the 
imaginary parts of higher-order diagrams describing the boson self-energies, resummed to all 
orders. This procedure has been used with great success in the past: indeed, the Z resonance 
can be described to very high numerical accuracy. However, in doing a Dyson summation of 
self-energy graphs, we are singling out only a very limited subset of all the possible higher-order 



'For differential distributions that do not involve an explicit phase, like e.g. dtot or dcr/(dcos 6*), there 
will be no interference between doubly- and singly- resonant diagrams. Consequently the non-doubly-resonant 
contributions are suppressed by T^/M^ sa 0.1% in those cases Q. 

™There are similar poles associated with diagrams containing internal Z propagators. These correspond to 
Z-pair production and are here considered as background to W-pair production. 
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diagrams. It is therefore not surprising that one often ends up with a result that retains some 
gauge dependence. 

For example it is very tempting to systematically replace — M^) by — M^ + iMF), 
also for < 0. Here F denotes the physical width of the particle with mass M and momentum 
q. This is the so-called 'fixed-width scheme'. As in general the resonant diagrams are not 
gauge-invariant by themselves, this substitution will again destroy gauge invariance. Moreover, 
the 'fixed-width scheme' has no physical motivation. In perturbation theory the propagator for 
space-like momenta does not develop an imaginary part. Consequently, unitarity is violated 
in this scheme. To improve on the latter the constant width could be replaced by a running 
one. This can, however, not cure the gauge non-invariance problem. At this point one might 
ask oneself the legitimate question whether the gauge breaking occurring in the 'fixed-width 
scheme' is numerically relevant or, like the gauge breaking in the LEPl analyses, negligible for 
all practical purposes. Of course, such a statement can only be made on the basis of an explicit 
comparison with truly gauge-invariant schemes in the full LEP2 energy range. We will come 
back to that point later on, after having defined a scheme that is gauge-invariant and reliable 
at LEP2. 

Below we will list a few ways to come to a gauge-invariant result and discuss their validity. 



One way to sidestep the gauge non-invariance problem is by simply multiplying the full 
matrix element by [q^ — M^]/[g^ — + zMF(g^)], which is evidently gauge-invariant [^, ^ . 
In this way the pole at = is softened into a resonance, at the expense of mistreating 
the non-resonant parts. It should be noted that when the doubly-resonant diagrams are not 
dominant, like at energies at and below the W-pair production threshold, this so-called 'fudge- 
factor scheme' can lead to large deviations 



The second possibility is the so-called 'pole scheme' |23|-|2^. In this scheme one decom- 
poses the complete amplitude, consisting of contributions from doubly-resonant diagrams 
(corresponding in lowest order to A^Born), singly-resonant diagrams R-, and non-resonant 
diagrams A^, according to their poles as follows: 
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- Ml, 



+ N{ki,ki,e) 



(5^ 



Here 6 stands generically for all angular variables which should be defined in such a way that 
their integration boundaries are independent of A;^ and k"^. Otherwise these angular variables 
would introduce an additional dependence on k'^ and k"^, and the correct pole terms could 
only be extracted after the angular integrations had been performed. This would complicate 
the pole decomposition and, in particular, would not be suited for a Monte-Carlo generator. 
Appropriate variables are e.g. the angles in the e+e~-CM system but not the Mandelstam 
variables. In equation (|58|) Ai is decomposed into gauge-invariant subsets originating from 
double-pole terms, single-pole terms, and non-pole terms (with repect to Mw). Introducing 
now the finite width only in the pole factors and not in the finite constant residues in brackets 
does not destroy gauge invariance. We note that different ways of introducing the finite width, 
e.g. constant or running, differ only by terms that are of higher-order and/or that are not of 
the double-pole type^. The same holds for different choices of the angular variables 6 (see 
section ^.4.41 ). 



A drawback of the 'pole scheme' is the fact that it is not defined below the W-pair production 
threshold and that it yields unreliable results just above this threshold (see section |3.4.4|) . 

Apart from yielding gauge-invariant results, the 'pole-scheme' decomposition also consti- 
tutes a systematic expansion according to the degree of resonance, i.e. in powers of rw/(Mw/3). 
Here the enhancement factor represents the infiuence of the nearby threshold on the expan- 
sion. Sufficiently far above the W-pair threshold and after imposing appropriate angular and 
invariant mass cuts, in order to reduce Z-pair and t-channel photon-exchange backgrounds, the 
cross-section for off-shell W-pair production is dominated by (or may even by defined by) the 
double-pole terms -R+_(M^, M^, 6*). At least at lowest order these may be related to on-shell 
W-pair production in the following way: 

R+-{Ml, Ml, e)= Y: -MeV-"-.w+w- X -^w+^/./2 X ■^w-./3/4' (59) 

where A^e+e-^w+w-; denote the matrix elements for the production of 

two on-shell W bosons and their subsequent decay into fermion-antifermion pairs. As such the 
'pole scheme' is a natural starting point for the systematic evaluation of higher-order corrections 
d, m (see also section p.4.4|) . 



As a third method, one may determine the minimal set of Feynman diagrams that is neces- 
sary to compensate for the gauge violation caused by the self-energy graphs, and try to include 



these [48, Hy]. This is obviously the theoretically most satisfying solution, but it may cause 



"There is no unique prescription for the propagator including the finite width. Instead of the constant term 
Mwrw one can also use a width depending on the invariant mass fc^ . This corresponds to different definitions 
of the (renormalizcd) W mass, which is fixed by the pole of the propagator. A popular choice is fc^Fw/Mw, 
which amounts to a shift in Mw by r^/(2Mw) ~ 26MeV relative to the propagator with the constant-width 
term. 
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an increase in the complexity of the matrix elements and a consequent slowing down of the 
numerical calculations. For the vector bosons, the lowest-order widths are given by the imagi- 
nary parts of the fermion loops in the one-loop self-energies. It is therefore natural to include 



the other possible one-particle-irreducible fermionic one-loop corrections 0, For the 
process e"'"e~ — > 4/ this amounts to adding the fermionic triple-gauge-boson vertex corrections. 
The complete set of fermionic contributions form a gauge- independent subset and obey all Ward 
identities exactly, even with resummed propagators . This implies that the high-energy and 
collinear limits are properly behaved. In contrast to all other schemes mentioned above, the 
'fermion-loop scheme' recommended here does not modify the theory by hand but selects an 
appropriate set of higher-order contributions to restore gauge invariance. To solve the problem 
of gauge invariance related to the width, we in fact only have to consider the imaginary parts 
of these fermionic contribution^. 

The 'fermion-loop scheme' should work properly for all tree-level calculations involving 
resonant W bosons and Z bosons or other particles decaying exclusively into fermions. This 
also includes, for instance, the hard-photon process e^e~ — >■ 4/+7, which requires in addition to 
fermionic vertex corrections also fermionic box corrections. For resonating particles decaying 
also into bosons, such as the top quark, or for calculating RC's to e"'"e~ — > 4/, which also 
involves bosonic corrections, the 'fermion-loop scheme' is not really suited. 

Although the latter scheme is well-justified in standard perturbation theory, it should be 
stressed that any working scheme is arbitrary to a greater or lesser extent: since the Dyson 
summation must necessarily be taken to all orders of perturbation theory, and we are not able to 
compute the complete set of all Feynman diagrams to all orders, the various schemes differ even 
if they lead to formally gauge-invariant results. In |Q another technique, the so-called 'pinch 



technique', has been introduced in order to construct a gauge-parameter-independent Dyson 
summation. Even if the 'pinch technique' yields gauge-independent results, it still contains 
some arbitrariness in the sense that one still has the freedom to shift gauge- independent parts 
that fulfill the Ward identities from the vertex corrections to the self-energies. For instance, it 



has been demonstrated in ||5J] that the 'background-field method' can be used to construct an 
infinite variety of such shifts, all representing (theoretically) equally well-justified schemes for 
resumming self-energies. 

Now it is a numerical question how much the predictions of different schemes differ. In 



5T| a detailed study has been given for the process e+e~ e~t'eud, a process that is highly 
sensitive to U{1) electromagnetic gauge violation. In this process the electron may emit a 
virtual photon, whose fc^ can be as small as m^: with a total centre-of-mass energy of ^/s 
available, we have a mass ratio of s/m^ = O{10^^), large enough to amplify even a tiny gauge 
violation in a disastrous wayQ. In table |^ we give the cross-section corresponding to the t-channel 
photon-exchange diagrams, responsible for the amplification of the gauge-breaking terms in the 
collinear limit. The results are given for two values of the minimum electron scattering angle 

°As the Ward identities are linear, we can separate the real and imaginary parts. 
^'This was noted already in [B5| , and investigated further in |M . 
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^min; displaying the effect of cutting away the dangerous colhnear hmit. It is clear that a naive 
introduction of a running width without a proper inclusion of fermionic corrections to the three- 
vector-boson vertex, which breaks U{1) electromagnetic gauge invariance, leads to completely 
unreliable results. The above-described gauge-invariant methods as well as the f/(l)-preserving 
'fixed- width scheme' numerically deviate by much less than Fw/Mw- Hence, a naive running 
width is not suited for LEP2, whereas a constant width, although SU{2) x f/(l) gauge breaking, 
might constitute a workable approach. As in the colhnear limit fc^ — the gauge-breaking terms 
originating from a naive running width are proportional to the dominant lowest-order graphs, it 
is possible to multiply the 7WW Yang-Mills vertex with a simple factor to successfully restore 
the U{1) gauge invariance. This factor is, however, certainly not universal. It will depend on 
the way the running width is introduced and on the process under investigation. Moreover, 
such a simple factor breaks unitarity and at high energies the full expression from the fermion 
loops is required for having a proper energy dependence. 



Scheme 






a no 


^min 10 


Fixed width 

Running width, no correction 
Fudge factor, with running width 
Pole scheme, with running width 
fermion-loop scheme 


.08887(8) 
60738(176) 
.08892(8) 
.08921(8) 
.08896(8) 


.01660(3) 
.01713(3) 
.01671(3) 
.01666(3) 
.01661(3) 



Table 9: Cross-section in different schemes for the t-channel photon-exchange diagrams of 
e"'"e~ —>■ e~z/eud. All schemes were computed using the same sample, so the differences are 
much more significant than the integration error suggests. 



3.4 Radiative corrections 

So far no complete treatment for the 0{a) corrections to off-shell W-pair production is available. 
Essentially only the initial-state photonic corrections, the final-state Coulomb correction, and 
the full hard process e"'"e~ — 4/+7 have been treated so far0. These are discussed in the 
following. The leading weak effects are normally taken into account through dressed lowest- 
order matrix elements, using and a(s).p In addition we describe a general strategy for the 
calculation of corrections beyond the lowest order using the 'pole scheme'. 

*An evaluation of all rcsummed one-particle-irreducible fermionic 0{a) corrections, in the context of the 
'fermion-loop scheme', is in progress |52[ . 

''The uncertainty associated with different theoretical definitions of s^, i.e. = 1 — Af^/Af| or = 
Tra{AM^)/{V2G^ M^), are found to be below 0.1%. 
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3.4.1 Initial-state radiation 



Most of the published calculations for corrections to off-shell W-pair production that have been 
done so far cover only part of the photonic corrections, mainly because these are easily treatable 
and constitute a large part of the RC's. From the discussion of the on-shell process we know 
that initial-state radiation (ISR) yields large corrections originating from the leading collinear 
logarithms. As will be discussed in appendix ^ these can be easily obtained by applying the 
structure-function method. Just in the same way as in the on-shell case the corrections to the 
off-shell cross-section are calculated by convoluting the lowest-order cross-section (in a certain 
gauge-invariant scheme) with the appropriate structure functions. In appendix ^ a detailed 
analysis is given of the theoretical uncertainties associated with the leading-log procedure. 

In the present-day Monte Carlos different ways of implementing the leading logarithmic 



corrections have been adopted One method involves solving the evolution equations for 
the structure functions numerically using techniques known from parton-shower algorithms. In 
this way soft-photon exponentiation and resummation of the leading logarithms from multiple 
hard-photon emission are automatically taken into account. Photons are generated according 
to the matrix elements in the collinear limit, in this way allowing for a finite pt kick to the 
photon. The second method involves a fully inclusive treatment of the radiated photons by 
folding the improved lowest-order cross-section with leading-log structure functions. So, no 
explicit photons are generated. The most recent development is a kind of merger of an explicit 
e"'"e~ — *• 4/7 Monte Carlo folded with structure functions, allowing for a consistent definition 



of observable and unobservable photon radiation (see section |3.4.3|). 



Another approach is to include the complete ISR. To this end one has to define it in a 
way that preserves the U{1) electromagnetic gauge invariance. This is non-trivial because of 
the presence of the t-channel diagram which involves a non-conserved charge flow in the initial 
state. One technique to circumvent this problem is the so-called current-splitting technique ||56| , 
which amounts to splitting the electrically neutral neutrino in the t-channel diagram into two 
oppositely flowing leptons each with charge one. One of them is attributed to the initial state to 
build a continous flow of charge, the other is attributed to the final state to do the same there. 
The modified Ue propagator leads to additional real and virtual initial-state photonic diagrams 
which render the ISR gauge-invariant. In this way one obtains the usual universal s-channel 
ISR known from LEPl plus additional non-universal contributions arising from the t-channel 
diagrams. The latter are non-factorizing with respect to the lowest-order cross-section, but 
are screened by a factor /c^/c^/s^, which automatically guarantees a unitary behaviour at high 
energies. Moreover they turn out to be numerically small at LEP2 energies. The universal ISR 
contains all leading logarithms and can be supplemented by the known universal higher-order 



terms IS? 



It should be noted that both the leading-log and the current-sphtting method leave out 
non-leading photonic corrections, for instance associated with radiation of photons off the 
intermediate W bosons. For on-shell W bosons we give in table ^ a comparison of exact and 
leading- log evaluations of the quantity / dE^ (da/dE^), needed for the average energy loss 
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leading-log results 


[GeV] 


[GeV ■ pb] 


0^ = s 

T 

1j 


0^ = s 

T 1 

— 1 


T 1 
Ij — 1 


161.0 


0.1377 ± 0.0001 


0.1436 


0.1379 


0.1375 


165.0 


3.619 ± 0.002 


3.792 


3.642 


3.607 


170.0 


10.120 ± 0.006 


10.651 


10.232 


10.089 


175.0 


17.437 ± 0.011 


18.431 


17.708 


17.403 


184.0 


30.882 ± 0.029 


32.873 


31.589 


30.902 


190.0 


39.562 ± 0.037 


42.224 


40.578 


39.595 


205.0 


59.126 ± 0.076 


63.581 


61.117 


59.322 



Table 10: The quantity / dE^ E^^ {da / dE^) needed for the average energy loss. The exact 
on-shell result from hard-photon radiation is given as well as the corresponding leading-log 
approximation, using L = log(Q^/m^) or L — 1. 



It is clear that a sufficiently accurate leading-log evaluation should be based on L — 1 
rather than L. Furthermore, the often-used scale choice = s leads to deviations of about 
3% at 190 GeV, which would translate into an error of ~ 60MeV on (E^). The scale choice 
Ql = 4Mw-E/(l + P), however, agrees with the exact result at the 0.3% level for all LEP2 
energies, leading to errors below lOMeV in (E^). As there is no obvious reason why these 
non-leading terms should be smaller in the off-shell case, some care has to be taken with the 
choice of a suitable leading-log scale. 



3.4.2 The Coulomb singularity 

Another potentially large photonic correction, not associated with the initial state, is due to 
the Coulomb singularity in the threshold region, which has been discussed for off-shell W-pair 
production in |^-|60|. It originates from the IR limit and at one-loop it emerges from a single 



IR-singular scalar three-point function and a related IR-singular scalar four-point function, the 
IR-singular part of which is just a scaled version of the one contained in the three-point function 
(see figure pi). The gauge invariance of the coefficients of these two scalar functions allows us 



"The (lowest-order) average energy loss (E^) can be obtained by normalizing to the lowest-order cross-section 
of table y. Dressing this lowest-order cross-section by LL structure functions, running couplings etc., only makes 
sense when the energy-weighted cross-section appearing in the numerator is treated in the same way! 
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Figure 8: Diagrams that contribute to the Coulomb singularity. 



to include the width in a gauge-invariant way. Thus one obtains a gauge-invariant correction 
factor to the lowest-order cross-section resulting from the doubly-resonant diagrams: 



with 



C^Coul = CT 



A 



CC3 

Born 



1 arctan 

TT 



■|/3m + A|2-^2 



M 



(60) 



^ - 2s{kl + kl) + {kl - P )2, 



1 - 4MVs, = M^- zMwFw - ie, 



\kl-kl 



(61) 



and — 7r/2 < arctan ?/ < 7r/2. Here /3 is the average velocity of the W bosons in their centre-of- 
mass system. Equation ( pDD only refers to the Coulomb singularity; the finite remnants of the 
above mentioned three-point function, containing also the IR divergence, have been left out. 
From equation (^) we obtain the usual on-shell Coulomb singularity for stable W bosons and 
/3 7^ by first going on-shell, i.e. A = and = Re/3|^ = and subsequently taking the 
limit Fw — > 0. Note that in this limit the otherwise negligible ie in becomes relevant. 



In contrast to the on-shell case there are various effects present in equation ( |60D that effec- 
tively truncate the range of the Coulomb interaction. The presence of a Coulomb singularity 
requires that /? be at least of the same order of magnitude as A and |/9m|, which is bounded by 

> 



the W width according to \I3m\ ^ y'^w/M^j. This confirms the intuitive argument that 
the Coulomb singularity is modified substantially by finite-width effects if the characteristic 
time of the Coulomb interaction (tcoui ~ l/[/3^Mw]) is of the same order as or larger than 
the typical decay time {tr ~ 1/Fw) of the W bosons. While the effect of the finite W width 
is contained in /Sm, the off-shellness shows up in A and in the difference — which 
effectively involves {k"^ + k"^ — 2M^\/s. This is in agreement with the argument that the 
off-shellness of the virtual W bosons will affect the Coulomb singularity, if the typical times 

for which off-shell bosons with four- 
It should be noted that as a 



l/\k. 



momenta k± exist is of the same order as or smaller than tcoui 
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result of the convolution with the Breit-Wigner functions from the W resonances, the quantity 
A plays a negligible role in an adequate description of the Coulomb phenomenon |5^. Con- 



sequently the at first sight deviating formulae of [^] and |Q constitute equally well-justified 
representations of this phenomenon. 

In the case of the total CCS cross-section, the Coulomb singularity gives rise to a correction 
factor which reaches its maximal value of ~ 5.7% at the nominal threshold and drops smoothly 
below and above threshold. While it amounts to 2.4% at ^/s = 176 GeV it is only 1.8% at 
190 GeV. The corresponding effect on the reconstructed W mass, resulting from the pronounced 
energy dependence, is a shift at the level of 5-10 MeV. Because of the fact that the Coulomb 
singularity is screened by the finite width of the W bosons and the off-shell effects, higher-order 
Coulomb corrections are not important for off-shell W-pair production [^-[^. Moreover, 



bound states of the two W bosons do not have the time to form because of the finite-width 
effects. The typical time scale needed for the formation of a bound state tform ~ 1/(q;^Mw) ~ 
234GeV~^ is much larger than the typical decay time tr = l/Fw ~ 0.5 GeV^^. 



3.4.3 The hard-photon process 



As mentioned before, in W-pair production the distinction between initial- and final-state 
radiation is not unique, unlike in the case of neutral particles such as the Z boson. In the matrix 
element, the universal leading logarithmic parts are easily separable, but the non-universal 
finite terms do not split naturally. An accurate calculation of the photon spectrum, beyond the 
leading logarithmic approximation, thus has to take into account initial-state radiation, final- 
state radiation, radiation off the W bosons, and various interference effects. These calculations 
have been performed for the final state consisting of 4 fermions and one photon ||62 
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However, the forward emission of many photons, described well by the structure functions and 
parton-shower algorithms described in section |3.4.1| and in appendix 0, is not included in these 
calculations. Recently efforts have been made to combine the two approaches [pR |BB|. We give 
here the most salient points of these papers. 

First we discuss the one-photon matrix element. This contains all graphs with two resonant 
W propagators, including radiation off the W bosons and the four-vector-boson interaction. 
The resulting non-leading terms are negative and tend to decrease the energy lost in initial- 
state radiation. On top of that one can also include the non-doubly-resonant graphs. The 
universal ones that contribute for all channels are totally negligible just as in the 

non-radiative process, whereas the t-channel graphs show the expected coUinear log(l — cos6') 
behaviour. 

The next step is to combine this hard matrix element, which is needed for large angles, with 

which gives a good 



3e 



the resummed leading-log structure function |65] or parton shower 
description of (multiple-) photon radiation at small angles. This is done by using the exact 
matrix element (convoluted with initial-state radiation) outside a cone around the incoming 
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20 25 
E^ [GeV] 

Figure 9: Photon energy spectrum in the one-photon approximation (I7) and the same com- 
bined with resummed structure functions (full) for the reaction e"'"e~ /i+i/^r"//^. 



electrons and positron^. Inside this cone one subtracts the leading logarithmic contribution 
from the hard matrix element, adds the tree- level matrix element, and subsequently convolutes 
everjd;hing with initial-state radiation. Of course the structure function or parton shower has 
to be restricted to radiate within the cone or virtuality cut-off. In the former case the scale 
Q"^ = s (1— cos6'c)/2 is used, with 9c defining the cone. 

In order to compare this approach with the one-photon matrix element on the one hand, and 
a purely leading-log description on the other, we have to define some kinematical observables. 
The most relevant ones for the W-mass measurement are the observable and unobservable 
photon energies. We define these with respect to the ADLO/TH set of canonical cuts defined 
in A photon which passes all cuts is called 'observable', and one that is combined with one 
of the beams is 'unobservable'. Photons close to final-state particles are not counted either way. 
The unobservable photon energy is the quantity most relevant to the W-mass reconstruction 
without explicit photons, whereas the observable photon energy gives an indication how well 
the resummed leading-log and one-photon matrix elements describe large-angle radiation. 

In order to get sensible results for these photon energies we will have to include an estimate of 
the unresummed soft and virtual corrections, which are not yet fully known. These corrections 
are only needed inside the cone as a weight of the structure function. As such the corresponding 
uncertainty partially cancels out in the average photon energies (E^) = J dE^ (da/dE^) E^/a. 
In the actual analysis presented below we circumvent the problem by leaving out the virtual 
corrections and tuning the non-leading-log part of the soft corrections (through the IR regulator 



'Instead of a cone also a cut-ofF on the virtuality can be used to define the region of multiple-photon 

radiation 1661. 
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leptonic 


semi-leptonic 


hadronic 


o-Q [pb] 


0.712 


4.485 


7.024 


< [pb] 


0.594 


3.732 


5.845 


^o+bkg [pb] 


0.595 


3.736 


5.853 










full 


0.275 


0.249 


0.229 


17 


0.365 


0.328 


0.302 


LL 


0.286 


0.260 


0.232 


with backgr 


0.275 


0.250 


0.229 


(Ef ^) [GeV] 








full 


1.227 


0.995 


0.760 


I7 


1.437 


1.167 


0.899 


LL 


1.238 


1.028 


0.801 


with backgr 


1.226 


0.994 


0.759 


^^unobs^ [GeV] 








full 


0.676 


0.664 


0.645 


17 


0.878 


0.873 


0.872 


LL 


0.666 


0.667 


0.664 


with backgr 


0.676 


0.664 


0.644 



Table 11: Cross-section for observable photons and the energy lost to observable and unobserv- 
able photons at y/s = 175 GeV. 



mass Air) in such a way that it cancels in the total cross-section against the non-leading-log part 
of the hard-photon corrections. Consequently, the resulting total cross-section is simply given by 
the lowest-order cross-section convoluted with leading logarithmic structure functions. As the 
initial-final and final-final interference is expected to be of order a^^^/M^^ (see section ^.4.41 ), 
this is a reasonable approximation. 



The results for the structure-function algorithm of [^ are given in table [TT] for leptonic, 
semi-leptonic, and hadronic final states^. For the cone 6c = 10° was chosen, but the results 
do not depend strongly on this parameter {9c = 5° only shifts the values by about 1%). In 
table |ll] we first give the observable tree- level non-radiative doubly- resonant cross-section (ctq) , 
the same convoluted with leading-log structure functions (ctq''), and with in addition the uni- 
versal non-doubly- resonant (background) graphs (o'Q'^+bkg)- -^^^ next entries we consider the 

"It should be noted that only the universal background diagrams were taken into account. 
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full calculation described above, the exact one-photon matrix element (I7), and the leading-log 
result (LL). We also give the full result including the universal non-doubly-resonant diagrams 
(full-|-backgr). The statistical errors on the cross-sections are (9(0.1%), but the differences 
(LL — full) and ((res -|- backgr) — res) were computed directly in this form and have relative 
errors of a few per cent and a few tens of per cent, respectively, on these differences. The 
statistical errors on the average photon energies are slightly larger, 0.3-0.5%. Note that we 
introduced an upper cut-off on the photon energy to avoid the Z peak. This influences only 
the observable average energy. 

One can see that an appreciable fraction (around one quarter) of the events will be ac- 
companied by photons observable in the ADLO/TH set of cuts. This is due to the excellent 
forward coverage {0^ > 1°) and electromagnetic calorimeter {E^ > 0.1 GeV) assumed in the 
canonical cuts. Reducing the angle to 10° this fraction still is around 20%, half of which also 
has an > 1 GeV. Neither the I7 matrix element nor the leading-log approximation give a 
satisfactory description of the observable photons. By radiating only one photon one misses 
the additional convolution with the structure functions, which scales down the one-photon 
cross-section as a result of the large, negative soft-photon effects. On the other hand, the 
leading-log approximation misses the negative effects from initial-final state interference and 
the radiation off the intermediate W bosons. Finally, given that the un-exponentiated large- 
angle contribution of the cross-section still is 20% of the total cross-section, the spectrum in 
the region E^ < 1 GeV should not be trusted to 1%, even not in the full calculation. Most of 
the cross-section here is, however, associated with the final state and hence does not infiuence 
the W-mass measurement. 

The observable photon energy is dominated by final-state radiation and hence not very 
interesting. The unobservable energy spectrum is much more independent of the specific final 
state. It is not completely independent due to the possibility of observable jets in the beam pipe 
(there is no angular cut on jets in the canonical cuts). The initial-state radiation associated with 
these events is sometimes combined with the final state by the canonical cuts, thus lowering the 
average energy. As the radiation off jets is not modelled correctly anyway, a jet-angle cut will 
have to be imposed to exclude this contamination. One sees that analysing observable photons 
separately reduces the average energy loss, and hence the size of the theoretical corrections to 
be apphed to the fitted W-mass. For the unobservable radiation the leading-log approximation 
is, as expected, quite good. Like in the total cross-section, the contribution of the universal 
non-doubly-resonant graphs is negligible. 

We conclude that the one-photon matrix element does not describe observable photons or 
unobservable photons well, whereas a leading logarithmic description is not accurate enough for 
large-angle photons, as it does not include the negative contributions associated with interfer- 
ence terms and radiation off the W bosons. Furthermore, a separate analysis of the events with 
an observable photon (around one quarter of all events with the canonical cuts) reduces the 
average energy loss considerably, thus reducing the uncertainties coming from this theoretical 
input. 
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3.4.4 General approach to radiative corrections within the pole scheme 



Finally, we discuss the way in which a reasonable approximation to the full radiative corrections 
to the process e"'"e~ 4/ can be computed. We do not consider QCD corrections. The 
perturbative QCD corrections to the final state are well-understood, as the main part is confined 
to the decay of a single W, which should be similar to the hadronic Z decays studied at LEPl. 
Perturbative QCD interference effects between the decay products of different W bosons are 
suppressed by a factor a"^/ (Nq — I). The related non-perturbative interference effects associated 
with the hadronization are discussed in section 13.51 



An electroweak 0{a) calculation would consist of the traditional three parts: one-loop 
graphs, soft-photon bremsstrahlung, and hard-photon bremsstrahlung. Only the sum is IR 
finite. The separation between hard and soft radiation is, as usual, defined with respect to a 
photon-energy cut-off -^min- Because of the finite W width one either has to use £^min ^ to 
ensure the validity of the soft-photon approximation in the soft-bremsstrahlung integrals (i.e., 
neglecting the effect of the photon momentum on the phase space and on the non-IR- singular 
parts of the matrix element), or one has to take into account finite photon energies in the 
Breit-Wigner resonances. Close to threshold, moreover, the limited amount of available phase 
space demands -Emin ^ 

The soft-bremsstrahlung integrals factorize into a simple multiplicative factor and the Born 
amplitude, and can easily be added, either using a 'fixed-width scheme' or the 'fermion- 
loop scheme'. Hard-photon radiation (for an arbitrary cut-off) has already been addressed 



in section 3.4.3. 



What remains are the virtual corrections. However, the full one-loop calculation of the 
virtual diagrams appears daunting. For the most simple final state, e"'"e~ — > fi^u^ud, there are 
3579 Standard-Model Feynman diagrams for massless fermions. This increases to 7158 when 
one electron is included in the final state, and reaches 15948 for the most complicated final state 
e"'"e~feZ7e. However, not all contributions are equally important. For instance the numerical 
significance of the non-doubly-resonant (background) diagrams is small in the Born and hard- 
photon calculations, especially if one requires that the event resembles W-pair production. So, 
we can assume that the one-loop corrections to these diagrams are even smaller. The problem 
then amounts to achieve a clean separation of the radiative corrections to W-pair production 



One way to achieve this separation is the 'pole scheme' introduced in section |3.3| for tree- 
level matrix elements. There have been various attempts to define one-loop corrections in 
this scheme ^^-|2^, which differ considerably. Since sofar no actual calculation for W-pair 



production has been completed in such a scheme, we restrict ourselves to a few comments. 
The idea behind the 'pole scheme' is to include a minimal part of the higher-order corrections 
needed to generate a finite width. By making a systematic expansion both in the coupling 
parameter a and in the width F (oc a), one can identify gauge- invariant contributions of 
progressively smaller infiuence on the final result. We will roughly sketch the method for a 
single, neutral particle. For factorizable diagrams, i.e. diagrams that factorize into corrections 
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to the production, propagation, and decay of the unstable particles, the amplitude to all orders 
can be Dyson resummed as follows p^, |67| 



fact 9 9 / ^ I 9 9 

pz _ \^pZ _ ^z 



E 



+ 



p2-M^ 1-S'(M' 



2^ ' 



(62) 



where m denotes the real mass, and M the complex mass given by — — S(M^) = 
0. The corrections to production and decay are contained in the function W, and S(p^) is 
the one-particle-irreducible self-energy. The variable 9 in the vertex-correction function W 



stands for other variables, which have to be chosen sensibly (see section p.3| ). The first term of 
equation (|62D does not have a pole and the residue at the pole = of the last term is gauge- 
invariant. In it has been shown how to derive the single-pole residue iy(M^, — S'(M^)] 
and the non-pole terms. The former consists essentially of the on-shell amplitude plus some 
terms of 0{r). It should, however, be noted that this type of decomposition only works when 
the on-shell limit exists. So, for instance below the W-pair production threshold the 'pole- 
scheme' method makes no sense, and as a result of that the procedure is not reliable just above 
threshold. 

When one turns to the production of charged unstable particles, additional problems arise 
related to the infra-red divergences that occur in the limit F — 0. Again these problems can 
be overcome 
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So far the above discussion only took into account factorizable diagrams. There are also non- 
factorizable diagrams, e.g. diagrams in which a photon connects the initial and final state or the 
decay products of different unstable particles. Those non-factorizable diagrams can give rise to 
double-pole contributions when the virtual photon becomes soft. Combined with the related soft 
bremsstrahlung it has been shown that the non-factorizable double-pole contributions cancel 
up to order aT/m in the fully inclusive total cross-section 
distributions this is in general not the case 



69, U^. For sufficiently exclusive 



3.5 Reconnection effects 

Nearly half of all W pairs decay hadronically. In the LEP2 energy range the average space- 
time distance between the W+ and W~ decay vertices is smaller than 0.1 fm, i.e. less than a 
typical hadronic size of 1 fm. Therefore the fragmentation of the W"^ and the W~ may not be 
independent, and this could influence the W-mass reconstruction. Here a short summary of the 
problem is given in a historical order. More can be found in the report of the W-mass group 

i- 

The problem is related to two different physical effects: colour reconnections and Bose- 
Einstein effects. 
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The colour-reconnection problem was first pointed out in ||72[ using the string-model ap- 
proximation. Assuming that the W pair decays into two quark- ant iquark pairs qq and QQ, 
respectively, it can either fragment into two strings stretched between qq and QQ., or into two 
strings stretched between qQ and qQ. In the three-colour world the probability for the second 
configuration was taken to be 1/9 in a very crude approximation. The properties of these two 
configurations are very different. While the first one gives 'classical' event multiplicities and a 
flat rapidity distribution, the second one has a small multiplicity and the particles are grouped 



at large rapidity values |72 



The fragmentation of the initial qq pairs into hadrons is conventionally described in terms 
of a perturbative parton cascade followed, at a later stage, by a non-perturbative hadronization 
phase. In perturbative QCD (PQCD) the influence of one W fragmentation on the other one 



is called 'interconnection'. Here the timing is very important P, |T3|: the gluon-emission time 
of high-energy gluons, Tg ~ iZ-E'g, is shorter than the W lifetime, rw ~ l/Pw ~ 0.5GeV~^. 
Therefore the gluons with Eg ^ Fw are emitted independently from the qq and QQ systems, 
while gluons with Eg < may feel four colour charges (=}- PQCD interference) |]73[. The 
non-perturbative hadronization (at distances ~ Ifm) is NOT independent. The influence of 
PQCD interference is shown to be not very important |73|, |7i[. It is suppressed by a colour 
factor 1/{Nq — 1) = 1/8 as compared to the total rate of double primary gluon emission and, in 
addition, only gluons with Eg < Pw are concerned; thus the total suppression factor is about 
100 compared to a naive instantaneous reconnection scenario with all events reconnected. It 
is much more complicated to estimate the influence of the fragmentation effects. In the Lund 



string model analogies to two types of superconductors have been used [^. The reconnection 
probability is taken to be either proportional to the space-time volume over which the W"*" and 
W~ strings overlap (type I superconductor) or to the probability that vortex lines cross (type 
II superconductor). Using these models the fragmentation error was estimated to be 30MeV 
on the W-mass measurement. This error is increased to 40 MeV by adding the estimated error 
coming from the perturbative effects as well as from the interplay between perturbative and 
non-perturbative effects. Also an extended version of this approach has been investigated |^, 
where additionally the space-time evolution of the parton cascade, multiple reconnections, and 
a finite vortex-core radius have been taken into account. 

In the space-time picture is not used, except for the timing of the gluon emission. 
Rather it is assumed that reconnections reducing the total string length are preferred, as indi- 
cated by PQCD experience. During fragmentation, in addition to the standard string connec- 
tions between quarks and gluons emitted from each W, also any other possible interconnection 
(recoupling) between partons emitted from different W bosons is possible with a probability 
that may be very different from 1/9. The reduced string length leads to a reduction in mul- 
tiplicity in the central rapidity region. In this model the shift in the reconstructed W-mass 
varies between 6 and 60 MeV ('shortest-string' version) or between 13 and 130 MeV ('random- 
reconnection' version) if the recoupling probability varies between 10 and 100%. Methods to 



measure this probability using data from LEP2 have been suggested [g| |75|. A similar model 



has been introduced in the Ariadne program leading to similar mass shifts Wt 
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The preliminary estimate of the colour reconnection in Herwig, based on a reduction of 
the space-time extent of clusters, gives mass shifts of the same order as the ones mentioned 
above §. 

The dependence of the shift in the reconstructed W-mass on the parameters of different 
models has been studied by the LEP experiments in more details 0. It appears that uncer- 
tainties using experimental procedures are not smaller than the ones obtained above. The issue 
whether diagnostic signals for reconnection can be found is still open. 

Thus the colour-reconnection effect may add a big systematic error to the W-mass mea- 
surement in the four-quark channel. The interesting aspect is that different approaches give 
uncertainties of the same order. 

It may be worthwhile to note that similar (but not identical) 'unconventional' colour con- 
nections can also appear inside a single Z system, as is discussed in p. The presence/absence 
of a signal in Z decays at LEPl could give an indication of the expected effects in the W^W~ 
system. 

Bose-Einstein correlations have been observed experimentally as an enhancement in the two- 
particle correlation function for identical bosons. In the case of a W pair decaying hadronically 
a possibility exists of Bose-Einstein correlations between particles that come from different W 
bosons [|7^]. A test of this effect has been made in a model based on Pythia and Jetset | |77| . 
The model gives large reconstructed positive mass shifts of the order of 100 MeV, rising with 
increasing c.m.s. energy and with decreasing source radius. The real shift may be smaller due to 
various damping factors not included in (which is intended as a 'worst-case' scenario), but 



effects of the order of 50 MeV could be expected within a large class of possible Bose-Einstein 
models. 

In the coming years both colour (re) arrangement and Bose-Einstein effects can and should 
be studied both theoretically and experimentally. Also the large statistics collected already at 
LEPl may help to study both these effects. 



4 Concluding Remarks 

In this report the present status of the theoretical knowledge of W-pair production and the 
related process of four-fermion production is reviewed. 

In lowest-order complete evaluations for all four-fermion final states exist, taking into ac- 
count all Feynman diagrams. The theoretical problem how to incorporate an energy- dependent 
width in the W and Z propagators without spoiling gauge invariance has been solved for the 
lowest-order cross-section. The numerical relevance of imposing electromagnetic gauge invari- 
ance is explicitly demonstrated. 



48 



As to the radiative corrections, the state of the art depends on whether on- or off-shell 
W-pair production, or four-fermion production is considered. Complete 0{a) RC's only exist 
for on-shell W-pair production. The dominant RC's that can be taken into account properly 
for all three types of processes comprise leading-log ISR, the leading weak effects related to 
and a{s), and the Coulomb singularity. The 0{a) RC's for four-fermion production obviously 
do not exist, and can at best be calculated using approximative techniques. Only taking along 
the above dominant effects, the total and differential cross-sections at 175 and 190 GeV are 
estimated to be known at the 1-2% level (based on on-shell experience). For the total cross- 
section at 161 GeV an uncertainty of about 2% is expected. Part of the latter uncertainty is 
due to the possible variation of the Higgs mass. This Higgs-mass dependence is largest for light 
Higgs bosons and is most pronounced at threshold. 

Another question, relevant for the reconstruction of the W mass, is the accuracy of the 
average emitted photon energy (Ery). At 175 and 190 GeV a precision of 10-20 MeV seems 
feasible, provided a proper scale is used in the structure-function part of the calculation. 



A Various Methods of Calculating QED Corrections 
A.l The structure-function method 

As pointed out in the previous sections, the virtual and real corrections reveal the presence of 
large logarithmic QED effects of the form aL/n = [a/ir) log(Q^/mg) with ^ m^. They 
arise when photons or light fermions are radiated off in the direction of incoming or outgoing 
light particles |7^, |79| , provided the momentum of the latter is kept fixed (exclusive) |j79l . In a 
massless theory these large logarithms would show up as collinear divergences (like in QCD), 
but in the SM the masses of the particles act as natural cut-offs. The fact that they can 
nevertheless give rise to sizeable corrections is due to the difference in scale between the mass 
of the radiating particle and its energy. This is illustrated by considering the propagator P of 
a light particle after photon emission 



1 qo'>rn 

P 




{q — ky — '^Qoko [1 — (1 — m? /2ql) cos 6] 



In the limit m ^ this propagator gives rise to a pole at cos 9 = 1. For finite m it yields 
large logarithmic terms of the form log(go/''^^) when the photon momentum is integrated over 
(inclusive photon). A consequence of the direct relation between the large QED logarithms 
and collinear divergences is that they are controlled by renormalization group equations and 
that they are universal, i.e. they are process- independent. They can be calculated using the 



so-called structure-function method [16|, taken over from QCD. This procedure also allows the 



inclusion of soft-photon effects to all orders by means of exponentiation. 
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Figure 10: Initial-state collinear radiation in e+e — > W'^'W . The unlabeled external lines 
represent an arbitrary number of undetected (inclusive) particles. 



The structure-function method is based on the mass-factorization theorem (see figure |T0]) 



^,dW^w^^ ^ic^ (63) 
dti aui at2 du2 Jo x+ Jo xt dti dui dt2 du2 

which links the cross-section cre+o-^w+w-X; where X indicates an arbitrary number of un- 
detected particles, to the hard scattering cross-section o"jj^w+w-x' which is free of the large 
collinear logarithm^. The indices i,j represent all light-particle transitions allowed by QED 
(e.g. 7, e^, n^, ■ • •). Furthermore we have defined 

ti = (p_-A;_)2-M^, ui = {p+-k_f ~ M^, 

t2 = ip+ - k+f - M^, U2 = (p- - k+f - M^. (64) 

Analogously their hatted counterparts are given by the same expressions with p± replaced by 
p± = x±p±. It should be noted that (^) allows the implementation of various cuts on the 
energies and angles of the produced W bosons, and that it can be used to extract more com- 
monly used distributions [like da / {dcos 9+)] by supplying the appropriate Jacobians [^]. The 



^If the final state consists of exclusively treated light particles (e.g. e+e~ 4/), large collinear final-state 
QED logarithms appear. These large logs can be treated in a way similar to the initial-state ones, provided 
proper account is taken of the fact that now p — p/x. 
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structure functions Tij describe 'mass singular' initial-state collinear radiation and represent 
the probability of finding in the parent particle j at the scale a particle i with fraction x 
of longitudinal momentum. So {1 — x) denotes the fraction of energy carried away by collinear 
radiation. All large collinear logarithms are contained in the structure functions, which depend, 
just like the hard scattering cross-section, on the mass-factorization scale Q^. These structure 
functions can be decomposed into a part containing the collinear logarithms and a non-log part 
according to 



X] 



a 



n=l ^ " ^ 1 



n n 



(x) U 



m=l 



-\non-log , 
ij ' 



n=l 



a 



(65) 



Analogously the hard scattering cross-section can be decomposed into the scale-independent 
Born cross-section plus non-log higher-order terms 



da, 



Born 



+ E 

n=l 



(66) 



For convenience we have dropped the explicit dependence on x± in this decomposition. This 
dependence enters via the reduced momenta p±. 

The so-called leading-log (LL) approximation consists in only taking along the terms oc 
(oL/tt)", which automatically means that the hard scattering cross-section a should be iden- 
tical to the scale-independent Born cross-section. In addition this Born cross-section could be 
dressed by the leading corrections that are not of the LL type, e.g. running couplings or a 
Coulomb factor, in order to improve the approximation. As these LL contributions constitute 
the most important higher-order QED effects, it is on most occasions sufficient to add the 
0{a'^) LL contributions to the full 0{a) results and to exponentiate the soft-photon distribu- 
tion. In this way the resulting cross-section becomes scale-dependent. A further simplification 
can be achieved by realizing that in e"'"e~ collisions the bulk of the QED corrections frequently 
originates from pure photon radiation, especially when the soft-photon contributions are dom- 
inant"". 



After these simplifications the relevant structure function takes on the form pT| , |52 | 

2^ _ exp(-|7j;/?exp + f/^s) /?exn -.a 1 



r^,L'-p(x,Q2 

1 

~ 422! 
1 



433! 



/3| 
Pi 



'){a, X, Q 

1 + 3x2 
1 — X 

l + x) 



r(i + i^exp) 



^(l-a:)^-/2-i-^/3H(l + x) 



log(x) + 4(1 + x) log(l — x) + 5 + X 



6 Li2(x) + 12 log2(l - x) - 37r^ 



™For more details concerning polarized structure functions and QED corrections that are not of the pure- 
photon-radiation type we refer to M . 
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X 



3 

- (1 + 8x + 3x^) log(x) + 6 (x + 5) (1 - x) log(l - x) 



+ 12 (1 + x^) log(x) log(l - x) - ^ (1 + 7x^) log2(x) + i (39 - 24x - 15x^ 



}• (67) 



Here Li2(?/) is the dilogarithm, 'Je the Euler constant, and T{y) the Gamma function (not to 
be confused with the structure functions). In terms of this structure function the total W- 
pair production cross-section including exponentiated LL QED corrections can for example be 
written as 



crLL,exp(s,Q^) = I ^ dz(l){2a,z,Q'^)ao{zs), (61 



where aolzs) denotes the (possibly improved) Born cross-section at the reduced CM energy 
squared zs. 

Note that some non-leading terms can be incorporated, taking into account the fact that 
the residue of the soft-photon pole is proportional to L — 1 rather than L for the initial-state 
photon radiation. Using 

r] = —L and /3 = — (L - 1) (69) 

71 TT 

we can identify a couple of popular options for including non-leading terms: 
• /3exp = (3, (3s = f3n = vE'¥3 ( 'MIXED' choice) 



Pexp = Ps = Pu = V 1 84] ('ETA' choice); the original Gribov-Lipatov form factor |85 



• /Sexp = (3s = (3n = (3 mm ('BETA' choice). 

The differences between these options are in the coefficient in front of the exponentiated soft- 
photon term and in the use of ?7 or /3 in the hard parts. 

In table |T2| we quantify for the CCS process e+e~ — > W+ W~ 4/ the uncertainty associated 
with these three different structure functions. We give the total cross-section without cuts, the 
radiative energy loss (E^), and the invariant-mass loss (My). The radiative energy loss is 
defined as 

(E^) = ^ J (2-X+ -x_)Eda (70) 
and the radiative invariant-mass loss as 

(My) = - / (l-x+x_)Eda. (71) 
a J 



As can be seen from the comparison between the third and fourth column of table |T^, the 
effect on the total cross-section due to the different coefficient in front of the exponentiated 
soft-photon term in the MIXED and ETA structure functions is of the order of 0.3-0.4%. This 
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MIXED 


ETA 


BETA 


BETA 


[GeV] 




to 


to C(/32) 


to 0{(3^) 


to 0{p^) 


175 


a 


13.2218(11) 


13.1772(12) 


13.1828(13) 


13.1790(11) 






1.1120(3) 


1.1115(3) 


1.1147(3) 


1.1146(3) 






1.1093(3) 


1.1090(3) 


1.1120(3) 


1.1119(3) 


190 


a 


16.3412(8) 


16.2829(12) 


16.2969(13) 


16.2939(8) 






2.1303(4) 


2.1300(5) 


2.1395(5) 


2.1391(4) 






2.1220(4) 


2.1216(5) 


2.1314(5) 


2.1314(4) 


205 


a 


17.1302(17) 


17.0696(12) 


17.0897(17) 


17.0869(13) 




{E,) 


3.1823(9) 


3.1817(9) 


3.2011(9) 


3.2019(7) 






3.1666(10) 


3.1655(9) 


3.1845(10) 


3.1854(9) 



Table 12: Effects of different structure functions on a (inpb), {E^) (in GeV), and (M^) (in GeV), 
for tfie CCS process e+e" W+W" 4/. 

different choice in the coefficient does not affect (E^) and (M^), since it is a factorized soft- 
photon contribution which largely cancels out in the ratios. By comparing the fourth and fifth 
column of table 0, one can estimate the amount of the so-called r] ^ (3 effect in the hard 
part of the structure function. The impact on the cross-section is small, i.e. of the order of 
0.1% (depending on the energy). The effect on {E^) and (M^) is 3, 10, 19MeV (at ^/s = 175, 
190, 205 GeV, respectively) and therefore not significant in view of the expected experimental 
accuracy at LEP2. From the last two columns of table 0, finally, we can infer that the effects 
of the 0{l3^) hard part of the structure function is completely negligible. 

At this point we would like to note that the leading effects related to the production of un- 
detected e"'"e~ pairs from conversion of a virtual photon emitted from the initial state (+ corre- 
sponding loop corrections) can be accounted for by replacing a in 0(a, x, Q^) by a[l + aL/ {6n)]. 
This takes into account the QED /3-function contributions appearing in the renormalization 
group equations. Also other fermion pairs can be included, but care has to be taken with the 
fact that the real pair-production process requires the virtual photon to be sufficiently hard. 

Another source of uncertainties is related to the fact that the scale is a free parame- 
terQ. As mentioned already in section p.3.3| , all 'natural' scale choices are roughly equal close 



to threshold. Hence, when using the structure-function method to calculate higher-order cor- 
rections [beyond C(a)] the dependence is negligible. This is, of course, not the case for LL 
0{a) corrections, as they are larger. So, in situations where the exact 0{a) corrections are not 
known and one has to resort to a LL approximation instead, as is the case for e"'"e~ — > 4/, the 
scale dependence is larger. 

^For a discussion of 'appropriate' scale choices we refer to 0. 
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A. 2 The parton-shower method 



The basic assumption of the QED parton-shower method is that the structure function of an 
electron (or positron) obeys the Altarelh-Parisi equation which can be expressed by the 
integral equation 



ree(x, Q') = nee(g^ QI) T^X, Q]) + £ ^ ^ee(g^ K^) ^ Pee(l/) T,,{x/y, K^) 

' ^ (72) 

in the leading-log approximation [Q. Here e is a small quantity specified later and Pcc{x) = 



1 + X — x). The function Ilee, which is nothing but the Sudakov factor, is given by 



nee(Q',g") =exp 



a 



2tt Jq'^ Jo 



dxPppix) 



(73) 



and denotes the probability that the electron (or positron) evolves from Q'^ to without 
emitting a hard photon. The scale is, like in appendix [A.l| , a free parameter (of order ml). 



Often it is chosen such that, when the dependence is integrated out, the factor j3 emerges 
rather than t]. 

Equation (|72|) can be solved by iterating the right-hand side in a successive way. Hence it is 
apparent that the emission of n photons corresponds to the n-th iteration. As such it is possible 
to regard the process as a stochastic one, suggesting the following algorithm for the photon 
shower ||8^: (a) Set Xb = 1. The variable Xb is the fraction of the light-cone momentum of the 
virtual electron (or positron) that annihilates, (b) Choose a random number ^. If it is smaller 
than Ilee{Q^,Ql), then the evolution stops. If not, one proceeds by finding the virtuality 
that satisfies ^ = Ilee{K^, Q^)- With this virtuality a branching takes place, (c) Fix x according 
to the probability Pee{x) between and 1 — e. Then Xb is replaced by XbX. Subsequently one 
should go to (b), replacing by K^, and repeat until the evolution stops. 

Once an exclusive process is fixed by this algorithm, each branching of a photon in the 
process is dealt with as a true process, that is, an electron with x, decays as e~(x, —K^) 
e~{xy, —K'^) +7(x[l — y], Qq). Here Qq is a cut-off to avoid the IR divergence. It is unphysical, 
so any physical observable should not depend on it. The momentum conservation at the 
branching gives —K^ = —K''^/y + Ql/{1 — y) + kf'/(|/[l — y]), which in turn determines k."^ 
from y, K^, K'^ . Hence the distribution can be taken into account in the simulation as well 
as the shape of x. This feature represents the essential difference between the structure-function 
method, which treats the photons inclusively, and the QED parton-shower method. Keeping 
in mind that y < 1 — e, the kinematical boundary y{K^ + Qq/[^ — y]) < K'^, equivalent to 

> 0, fixes e to e = Ql/K''^, since strong ordering {K'^ -C K''^) is expected. 

The above description of the algorithm represents the "single-cascade scheme" . This implies 
that only the e~ or the e^ is able to radiate photons when the axial gauge vector is chosen along 
the momentum of the other initial-state particle, namely e^ or e~. For programming purposes, 
however, it is convenient to use a symmetrization procedure (the so-called "double-cascade 
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scheme") to ensure the symmetry of the radiation with respect to the electron and positron 



90| , |91| . The so-obtained QED parton shower can then be combined with the matrix element 



of any hard process initiated by e+e annihilation. 



The Altarelli-Parisi equations can also be solved by a different algorithm which yields 
equivalent results for electromagnetic radiative corrections. 

Starting from the observation that the most important corrections come from multiple soft- 
photon emission, the photon-number distribution can be approximated by a Poisson distribution 
with average photon number 

n^ = ^\og(9l] ['~'dxPUx). (74) 

The technical infra-red cut-off e will drop out of the physical sum over soft photons. After 
the number of photons has been chosen from the Poisson distribution, their transverse and 
longitudinal momenta can be generated properly, ordered according to the l/{p± ■ k) poles and 
splitting functions, respectively. 

Using momentum conservation at each branching, the correct energy loss and pt boost 
for the hard scattering is obtained. This algorithm has been shown to yield results that are 
consistent with the structure-function formalism and provides a realistic approximation for the 
photonic distributions. 



A. 3 YFS exponentiation 

Let us first answer the basic question: what the exclusive Yennie-Frautschi-Suura (YFS) expo- 
nentiation is and what it is notl YFS exponentiation is, in one word, a technique of summing 
up all IR singularities for a given arbitrary process to infinite order. The formal proof in ||68[ is 
based on Feynman- diagram techniques and goes 'order by order'. The YFS technique is not by 
any means bound to the leading-log summation technique and/or approximations. It is appli- 
cable to arbitrary stable particles (with arbitrary mass and spin) in the initial and final state. 
The YFS summation is inherently exclusive, i.e. all subtractions/summations of IR-singular 
contributions are done before any phase-space integrations over virtual- or real-photon four- 
momenta are performed. This means that the Monte Carlo technique can be used to integrate 
over the multiple real-photon phase space. The first practical solution of this kind was given in 



93| . Here the QED matrix element that enters the YFS exponentiation consists of two parts. 
The 0{a) part is taken to be exact, i.e. from Feynman diagrams, whereas a LL approxima- 
tion is used to write down an economic ansatz for the 0{a^) part. The YFS-exponentiation 
procedure, which involves a subtraction of the IR part of the matrix element, knows nothing 
about the origin of the matrix element (ansatz or Feynman rules), and it will fail if the ansatz 
has the wrong IR limit. The above subtraction procedure is done on the fully differential dis- 
tribution before phase-space integration. Hence it should be clear that a LL ansatz for the 
matrix element in which the emitted photons have zero transverse momenta is not possible. 
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Contrary to the LL techniques described in appendices A.l and A.2| , in the YFS approach the 
parameter (3 = 2 {a / n)[\og{s / ml) — 1] results from the integration over phase space and there is 
no freedom to adjust it. There is no discussion "do we have —1 in the definition of the leading 
logarithm or not". The —1 is mandatory, because otherwise the IR singularities do not cancel. 
The orthodox YFS exponentiation technique |^ is essentially 'order by order' and as such it 
does not sum up all the LL corrections to infinite order, something which is in principle possible 
for the LL techniques described in appendices |A.1| and |A.2| . 



As the exact 0{a) matrix element for e^e — > 4/ is not yet known, the present implemen- 
tation of YFS exponentiation for this process relies on a pure LL ansatz |^4|. But this will, of 



course, change as soon as an adequate approximative 0{a) calculation becomes available. 
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